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INTRODUCTION 


Terms  of  Reference 

This  Supplement  to  the  report  on  Program  for  Reclamation  of  Surface 
Aquifer,  Rocky  Mountain  Arsenal,  has  been  prepared  by  the  firm  of  C.  E  o 
Jacob,  Groundwater  Consultants,  Northridge,  California,  acting  as 
Architect-Engineer  under  Contract  No.  DA-25-066-eng-636l . 

Under  this  contract  the  Architect-Engineer  is  authorized  to  perform 
the  following  services; 

(Article  l-B  -  Statement  of  Services) 

"1.  Advise  this  office  on  collection  of  water-level  and 
water-quality  data  in  the  field  and  also  on  the  evaluation 
and  study  of  this  data. 

"2.  Initiate  mathematical  analyses,  model  studies  or 
numerical  calculations  of  the  following  problems: 

a.  Dispersion  of  salt  water  and  other  contaminants 
reaching  the  water  table  from  infiltration  or  evaporation  ponds. 

b.  Hydraulics  of  thin  unconfined  groundwater  systems 
on  curved  and  sloping  bottomSc 

c.  Influence  of  anisotropy,  inhomogeneties  and  strati¬ 
fication  upon  flow  of  groundwater  in  such  aquifers  and  upon 
dispersion  of  contaminants  therein. 

d.  Hydraulics  of  line  of  reclamation  wells  arranged  down- 
gradient,  or  cross-gradient. 

e.  Hydraulics  of  other  arrays  of  reclamation  wells. 

f.  Hydraulics  of  other  feasible  reclamation  schemes 
such  as  trenches. 

g.  Dispersion  of  distributed  contaminant  upon  attempted 
withdrawal  through  various  reclamation-well  or  trench  systems, 
also  residual  of  contamination  under  each  system  and  distribution 
thereof. 

"3.  Test  Drilling.  Advise  on  most  efficient  spacing  of  test 
borings  to  define  bedrock  and  outline  extent  and  properties  of 
aquifer.  The  positioning  of  these  borings  will  be  done  where  possi¬ 
ble  to  permit  their  use  in  pumping  tests  of  existing  or  specially  con- 


constructed  pumping  wells, 

"4,  Geophysieol  exploration.  Advise  on  interpretation  of  geo¬ 
physical  surveys  conducted  to  map  bedrock  between  boringS/  whether 
by  electrical  resistivity  or  by  seismic  methods. 

"5,  Piimpino-test  methods.  Outline  procedure  and  techniques  for 
determining  character  of  aquifer  by  well  pumping  tests.  Insofar  as 
possible,  tests  will  be  run  first  using  the  existing  production  wells, 
such  as  abandoned  farm  wells  still  equipped  with  pumps  in  nearby 
areas.  Design  of  pumping-test  layouts  requiring  specially  constructed 
production  and  observation  wells  will  be  planned  by  Architect-Engineer 
for  the  Omaha  District. 

•'6.  Leakage  of  ponds.  Advise  on  methods  of  testing  Pond  F  and 
other  ponds,  if  desired,  for  water-tightness.  This  may  be  done  by 
restricting  inflow  and  observing  water-level  decline,  or  by  using 
tracers  if  possible. 

"7.  Water  balance.  Undertake  to  reconcile  best  estimates  of 
historic  and  present  inflow  to  groundwater  system  and  displacement 
of  contaminated  water. 

"8.  -*^alt  balance.  On  the  basis  of  latest  geologic  data  and 
groundwater  data  plus  distribution  of  salts  revealed  by  program  of 
sampling,  make  an  inventory  of  various  signlgicant  ions  (such  as 
chloride).  These  salt  balances  will  be  reconciled,  if  possible,  to 
the  aforementioned  water  balance.  Attempt  will  be  made  to  explain 
the  past  and  present  contamination  in  the  light  of  the  known  character¬ 
istics  of  the  aquifer  and  of  its  geologic  structure  as  revealed  by  test 
drilling  and  sampling. 

"9.  Restricting  the  flow  at  arsenal  boundary.  Undertake  studies 
of  various  methods  proposed  to  restrict  flow  of  contaminants  under¬ 
ground  at  Rocky  Mountain  Arsenal.  These  analytical  studies  may  be 
supplemented  by  model  studies  or  by  electronic  analog  computations. 

"10.  Reclamatinn  nf  nguifer.  Different  schemes  for  removing  con. 
taminated  groundwater  through  wells  or  open  excavations  will  be 
evaluated.  The  possibility  of  injecting  fresher  water  after  pumping 
the  contaminated  water  to  accomplish  dilution  of  the  water  in  the 
ground  will  be  studied.  This  will  require  analytical  and  electronic 
analog  studies  of  the  further  dispersion  of  the  contaminants  upon 
attempted  withdrawal." 


Conferences 


Conferences  between  the  Architect-Engineer  and  personnel  of  the 
Omaha  District  and  of  the  Missouri  River  Division  of  the  Corps  of  Engineers 
were  held  in  Omaha  on  the  following  dates:  6  July,  9  August,  20-21  Octo¬ 
ber,  and  6  December,  1960.  During  these  conferences  the  Architect-Engineer 
advised  the  Corps  of  Engineers  on  Items  1 ,  3,  4,  5  and  6  under  Article  l-B  of 
the  contract. 

The  initiation  of  mathematical  analyses  and  numerical  calculations, 
under  Item  2,  is  covered  in  Appendices  E  through  I  in  this  Supplement. 

The  outlining  of  the  procedures  and  techniques  of  pumping  tests,  under 
Item  5,  is  covered  in  Appendices  B  and  C  of  this  Supplement. 

The  water  balance  and  salt  balance,  under  Items  7  and  8,  were  dis¬ 
cussed  in  the  Preliminary  Report  and  are  reviewed  in  this  Supplement. 

Methods  of  restricting  flow  of  contaminants  at  the  arsenal. bOUfKlary  ond 
of  reclaiming  the  aquifer,  under  Items  9  and  10,  are  covered  in  the  five  alter¬ 
native  schemes  outlined  in  the  report  itself.  These  schemes  were  considered  in 
detail  during  the  conference  of  6  December  1960  in  Omaha. 

Preliminary  and  Progress  Reports 

A  Preliminary  Report  on  Ground  Water  Problem,  Rocky  Mountain  Arsenal, 
was  written  1  August  1960.  That  report  contained  one  table  and  one  plate  and 
13  pages  of  text,  reproduced  in  blue-line  prints.  Progress  reports  in  the  form 
of  Memoranda  of  Conferences  were  prepared  9  August  1960  and  21  October  1960, 
and  were  reporduced  in  blue-line  prints. 


AQUIFER  CHARACTERISTICS 


The  hydraulic  characteristics  of  the  aquifer  were  determined  by  a 
number  of  well  tests  run  In  Area  A  (the  contaminated  area  extending  north¬ 
ward  from  the  northwest  boundary  of  the  Arsenal).  Table  S-1  gives  a 
summary  of  these  formation  constants.  In  the  first  column  the  kind  of  test 
is  Indicated  and  in  the  second  the  number  or  name,  of  the  wells  tested.  In 
column  3  values  of  hydraulic  conductivity  (K)  in  feet  per  day  are  given  and 
in  column  4  values  of  the  storativity  (S).  Pertinent  remarks  are  given  in 
column  5. 

For  the  locations  of  the  wells  tested  see  Figures  D-la  and  D-lb  in 
Appendix  D. 

Ten  wells  shown  on  Figure  D-la  were  given  drawdown-recovery  tests. 

The  results  were  doubtful  in  two  cases,  the  Marty  Well  and  the  Matsumoto 
Well  2.  The  hydraulic  conductivity  determined  by  the  remaining  eight  tests 
ranged  from  a  low  of  1 ,040  to  a  high  of  2,640  ft/day.  The  average  of  the 
eight  tests  was  1 ,483  ft/day.  (We  shall  round  this  off  at  1 ,500  ft/day  for  the 
purpose  of  calculation.) 

The  data  on  the  second  line  of  Table  5-1  refer  to  the  interference  test 
of  the  Aden  V/ell.  For  further  details  see  Tables  D-4  and  D-5  in  Appendix  D. 

According  to  the  logs  of  the  six  piezometers  or  observation  wells  (see 
Figure  D-47b)  the  groundwater  is  generally  unconfined  in  the  vicinity  of  the 
Aden  Well.  The  logs  of  three  of  the  observation  wells  show  "gravelly  sand" 
or  "sandy  gravel"  at  the  level  of  the  water  table,  whereas  two  of  the  other 
three  show  "clayey  silt"  or  "clayey  sand".  In  the  sixth  observation  well  (P-6) 
about  five  or  six  feet  of  sandy  clay  is  logged  at  and  just  above  the  water  table. 
Despite  these  variations  shown  in  the  logs,  a  condition  of  unconfinement  is 
reflected  in  the  variable  and  high  values  of  apparent  storativity  obtained  from 
this  interference  test. 

From  the  average  drawdown  and  recoveries  observed  in  the  six  observation 
wells  the  hydraulic  conductivity  is  calculated  to  be  about  1,230  ft/day.  The 
initial  value  of  apparent  storativity  is  0.056  or  5.6  percent.  By  extrapolating 
the  trend  of  apparent  storativity  it  is  estimated  that  the  maximum  value  would 
be  on  the  order  of  0.35  or  35  percent. 

The  data  on  the  third  line  of  Table  S-1  relate  to  the  interference  test 
on  the  Marty  Well.  According  to  the  logs  plotted  on  Figure  D-48b,  at  the  time 


TABLE  S-1 


SUMMARY  OF  FORMATION  CONSTANTS  FROM  WELL  TESTS 


[1] 

[2] 

[3] 

[4] 

[5] 

Kind  of 
test 

Well(s) 

Hydraulic 

conductivity 

K/[ft/day] 

Stora- 

tivity 

S 

Remarks 

Drawdown- 

recovery 

eight 

1,500 

— 

Ten  tests  in  all; 
two  doubtful 

1 nterference 

Aden 

1  ,230 

0.056 

Avg.  initial  value 

.35 

Max.  ultimate  value 

1 nterf erence 

Marty 

1,330 

.085 

S  fairly  const,  in  time 

1 nterf erence 

Powers 

1 ,090 

.10 

S  fairly  const,  in  time 

of  the  test  the  water  table  stood  In  "sand"  or  "gravelly  sand"  in  all  six 
observation  wells. 

The  average  hydraulic  conductivity  from  the  drawdown  and  subsequent 
recovery  observed  in  all  six  observation  wells  is  1,330  ft/day.  (See  Table 
D-6.)  The  storativity  remained  fairly  constant  throughout  the  period  of  the 
test  and  averaged  about  .085  or  8.5  percent.  Undoubtedly  the  storativity 
would  increase  upon  prolonged  pumping  or  drainage. 

The  data  on  the  fourth  line  of  Table  S-1  relate  to  the  interference  test 
of  the  Powers  Well.  Logs  of  four  observation  wells  drilled  near  the  Powers 
Well  are  shown  graphically  on  Figure  D-23b.  In  three  of  those  four  wells 
"clay"  or  "silty  clay"  was  encountered  at  or  just  above  the  water  table  and 
in  the  fourth  one  a  layer  of  "silty  sand".  At  the  time  of  the  test  the  water  table 
stood  in  "silty  sand"  or  "gravelly  sand"  beneath  these  semi-impervious  layers. 

The  average  hydraulic  conductivity  was  1 ,090  ft/day  and  the  apparent 
storativity  about  0. 10  or  10  percent.  The  storativity  was  fairly  constant  in 
time  over  the  period  of  drawdown  and  also  over  the  subsequent  period  of  recovery. 

In  summary,  for  determining  scale  ratios  for  digital  or  analog  computation 
we  shall  use  1,500  ft/day  for  the  hydraulic  conductivity  and  0.20  or  20  percent 
for  the  average  storativity. 
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WELL  CHARACTERISTICS 

Table  D-1  summarized  the  drawdown-recovery  tests  run  in  ten  wells  near 
Rocky  Mountain  Arsenal,  The  overage  depth  of  the  ten  wells  Is  37,2  feet.  The 
most  common  casing  diameter  is  48  inches  and  the  most  common  material  concrete. 
The  eight  wells  whose  tests  were  successful  had  an  average  discharge  of  about 
215  gpm  or  41,400  fr  /day.  The  average  hydraulic  conductivity  (see  column  [12] 
of  Table  D-1)  is  1,483  ft/day  or  approximately  1,500  ft/day. 

As  a  result  of  the  step-drawdown  test  on  the  Powers  Well  !t  was  concluded 
that  "for  practical  purposes,  over  the  range  of  discharge  involved  here,  when 
dealing  with  wells  similar  in  design  to  the  Powers  Well,  the  well  loss  may  be  neg¬ 
lected."  (See  Page  D-1 1 .) 

For  the  purpose  of  determining  scale  ratios  for  analog  or  digital  computation, 
then,  we  may  use  the  effective  well  radius  determined  from  the  drawdown-recovery 
tests,  its  unnaturally  small  value  making  allowance  for  the  transitional  flow  near 
the  well. 

Referring  to  column  [15]  of  Table  D-1 ,  it  is  seen  that  the  mean  value  of  the 
logarithm  of  the  quantity  (S  r^  /h*)  in  feet  is  0.72-5.  Taking  for  the  value  of 
h*  an  overage  initial  depth  of  flow  of  10.3  feet  (the  average  for  the  eight 
successful  tests),  the  logarithmic  mean  value  of  the  product  (S*r  )  is  about 
5.5  X  10*^  ff^.  This  is  the  figure  that  we  ^all  use  in  determining  scale  ratios 
for  electronic  analog  or  for  other  models  of  the  wells. 


CHARACTERISTICS  OF  CUT-OFF  TRENCH 


Figure  S-1  shows  a  series  of  water-table  profiles  upstream  and  downstream 
from  a  cut-off  trench  excavated  to  the  bottom  of  an  inclined  aquifer.  These 
profiles  were  obtained  by  numerical  integration  of  Equation  G-3  in  Appendix  G. 
To  interpret  them  in  physical  terms  it  is  necessary  to  have  an  estimate  of  the 
scale  factors  involved. 

If  the  initial  depth  of  flow  is  about  10  ft  on  the  overage  and  if  the 
slope  of  the  water  table  is  about  1/260  (See  Plate  7),  then  the  physical  dis¬ 
tance  (horizontal)  is  related  to  the  corresponding  nondimensional  space  coor¬ 
dinate  as  follows:  x=§  10  ft/( 1/260)  =  ^  2,600  ft.  Similarly  the  head 
measured  above  the  inclined  bottom  of  the  aquifer  is  related  to  the  corresponding 
nondimensional  coordinate  as  follows:  (h  -  z)  =  10  ft  (^“f  ).  (See  Appendix  E 
and  Glossary  for  meaning  of  symbols. 

The  physical  time  is  related  to  the  corresponding  nondimensional  time 
coordinate  os  follows:  t  =7''(Sh^/oc^K)  =  (0.20  x  10  ft  x  67,600/1 ,500  ft/day). 
Hero  we  have  assumed  that  on  the  average  S  =  20  percent  and  K  =  1 ,500  ft/day. 
It  turns  out,  then,  that  t  =  (90  days)  T'. 

Referring  again  to  Figure  S-1,  it  is  seen  that  after  about  0.6  units  of  time 
on  the  T-scale  the  gradient  downstream  from  the  trench  is  reduced  to  zero. 
Thereafter  the  flow  from  downstream  back  into  the  trench  is  zero,  and  provided 
the  head  or  water  level  In  the  trench  is  kept  at  or  below  the  bottom  of  the  for¬ 
mation,  there  will  be  zero  flow  from  the  trench  into  the  aquifer  downstream. 


FIG.  S-1 


WATER  TABLE  PROFILES  NEAR  CUT-OFF  TRENCH 
IN  INCLINED  AQUIFER 
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WATER  BALANCE 

A  check  may  be  made  of  the  rate  of  movement  of  the  contaminated  ground- 
water  and  of  the  time  that  might  have  been  required  for  it  to  traverse  the  distance 
from  Reservoir  A  to  the  South  Platte  Rivero  Water-table  contours  on  Plate  2  and 
the  profile  on  Plate  7  show  an  average  water-table  slope  or  hydraulic  gradient  of 
about  20  ft/mi.  The  average  hydraulic  conductivity  from  the  drawdown  -  recovery 
tests  is  about  1 ,500  ft/day.  Thus,  assuming  the  overage  porosity  is  35  percent, 
the  average  velocity  of  movement  is  as  follows 

V  =  (  20  X  1 ,500  ft/day)/(0.35  x  5,280)  =  16  ft/day 

At  this  rate,  without  dispersion  it  would  take  the  water  about  320  days 
to  traverse  a  distance  of  one  mile.  To  traverse  the  6.4  miles  distance  from 
Reservoir  A  to  the  River  it  would  take  about  2,100  days  or  5.7  years.  Actually 
the  most  highly  contaminated  water  has  not  yet  reached  the  river.  (See  Plate  4.) 

The  foregoing  calculation  applies  to  water  of  average  chloride  concentration 
(about  1,100  ppm).  By  longitudinal  dispersion  of  the  advancing  front  contaminated 
water  of  lower  concentration  reaches  the  river  much  sooner.  However,  by  transverse 
dispersion  of  the  advancing  front  the  advance  is  delayed. 

Although  evidently  the  contaminated  water  has  been  migrating  over  a  period  of 
about  17  years,  a  more  detailed  water  balance  than  is  now  possible  would  probably 
show  that  much  of  the  contaminated  water  has  gone  into  storage  by  lateral  dispersion 
on  either  side  of  the  advancing  thread  or  stream  flowing  from  Reservoir  A  to  the 
River.  The  influence  of  farm  irrigation  wells  in  helping  to  spread  the  contamination 
laterally  must  also  be  considered. 


SALT  BALANCE 


Estimates  of  volumes  and  concentrations  of  contaminated  water  in  the 
aquifer  under  and  near  Rocky  Mountain  Arsenal  given  In  Table  S-1  may  be 


irized  as  follows: 

On  the  Arsenal 

436  million  cu  ft  averaging 

1  /  680  ppm  chloride 

Off  the  arsenal 

- 

in  Area  A 

301 

720 

In  Area  B 

245 

560 

Total  982  million  cu  ft  averaging  1,  ICX)  ppm  chloride 

We  may  estimate  the  quantities  of  inflow  of  given  concentrations  needed 
to  produce  the  foregoing  estimated  volumes  of  contaminated  watero  These 
symbols  will  be  used: 

V  =  volume  of  water 

'Y  =  specific  weight  of  water/  weight/volume 
c  =  concentration  of  ion 
t  =  elapsed  time  of  flow 
Q  =  rate  of  flow,  volume  per  unit  time 
W  =  weight  of  ion  in  water  volume  V 

The  relationship  between  the  weight  of  any  given  ion  in  a  given  volume 
of  water  is  given  by  W  =  ctY.  Dividing  both  sides  of  the  equation  through 
by  time  we  get  the  following  relationship  W/t  =  c'vV/t  =  c  ^Q.  Solving  for 
the  rate  of  flow  of  water/  or  volume  per  unit  time/  we  get  the  following 
Q  =  W/cirt. 

To  make  a  rough  check  on  the  chloride  balance  it  will  be  assumed  that 
water  has  been  entering  the  water  table  and  contaminating  the  groundwater  over 
a  period  of  17  yearS/  or  about  6/000  days.  The  specific  weight  of  water  will  be 
taken  to  be  62.4  Ib/ft^.  From  the  foregoing  table  the  volume  of  contaminated 
water  Is  seen  to  be  about  7.34  billion  gallons  or  about  1  billion  cubic  feet. 
Thus/  W  =  crV  =  1/100  X  10“^  X  62.4  Ib/ft^  x  .980  x  loV  =  6.7  x  10^  lb. 
With  an  average  chloride  concentration  of  1/100  ppm/  the  total  weight  of 
chloride  ion  Is  estimated  to  be  67/000/000  pounds. 

Assuming  that  the  average  chloride  content  of  the  contaminating  fluid 
were  4/000  parts  per  million/  the  average  Inflow  rate  to  the  water  table  needed 


to  produce  the  estimated  volume  of  contaminated  water  would  be  as  follows, 

(In  this  calculation  the  initial  chloride  content  of  the  groundwater  of  the  area^ 
thought  to  be  on  the  order  of  50  ppm,  is  neglected.) 

Q  =  W/c  r  t  =  (6.7  X  10^lb)/(4,000  x  10"^  x  62.4  Ib/ft^  x  6,000  day) 

=  4.5  X  10^  ft^/day  =  235  gpm  ± 

Thus  it  is  seen  that  the  average  input  to  the  water  table  of  water  having 
an  overage  chloride  content  of  4,000  ppm  would  have  had  to  be  about  235  gpm. 

If  the  chloride  content  had  been  only  2,000  ppm  on  the  average,  the  average 
rate  of  Inflow  would  have  had  to  be  about  470  gpm. 

Petri  [1956,  p.  48]  reports  as  follows:  "Although  liquid  wastes  have  been 
discharged  into  the  disposal  reservoirs  in  Sections  26,  35,  and  36  since  1942,  no 
records  are  available  of  the  amounts  of  liquid  discharged.  Measurements  made  in 
June  1955  by  the  Ralph  M.  Parsons  Co.  ****  indicate  that  at  that  time  the 
overage  input  at  the  disposal  reservoirs  was  about  570  gpm. 

Obviously  the  volume  rate  of  input  of  contaminated  water  to  the  water 
table  has  continuously  been  less  than  the  volume  rate  of  input  to  the  disposal 
reservoirs.  However,  as  relatively  little  chloride  is  stored  In  the  abandoned 
disposal  reservoirs  or  retained  in  the  subsoil  of  their  bottoms,  the  rate  of  input 
of  chloride  to  the  water  table  has  probably  equalled  on  the  average  the  rate  of 
input  of  chloride  to  the  disposal  reservoirs.  Offhand  it  would  appear  that  the 
concentrations  and  quantities  of  Inflow  in  the  past  have  probably  fallen  within 
the  limits  set  by  the  foregoing  calculation. 
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THEORY  AND  PROCEDURES  OF  WELL  TESTING 
Definitions 

In  this  report  the  standards  of  the  Soil  Science  Society  of  America  and 
of  the  American  Petroleum  Institute  are  accepted,  and  the  following  definitions 
of  technical  terms  are  adhered  to. 

Permeability.  The  property  of  a  porous  water-bearing  material  that  relates 
to  its  ability  to  conduct  water  under  a  hydraulic  gradient  is  known  as  its  permea¬ 
bility.  The  author  uses  the  definition  of  this  term  accepted  in  the  pertoleum  in¬ 
dustry.  Thus,  the  permeability  is  the  horizontal  volume-flow  per  unit  area  per 
unit  time  per  unit  horizontal  pressure  gradient  of  a  fluid  of  a  unit  viscosity  and 
unit  specific  weight.  The  flow  may  be  expressed  in  cubic  centimeters  per  square 
centimeter  per  second,  the  horizontal  pressure  gradient  in  atmospheres  per  centi¬ 
meter,  the  viscosity  in  centipoises,  and  the  specific  weight  in  grams  per  cubic 
centimeter.  With  this  choice  of  units,  the  unit  of  permeability  is  called  the 
"darcy''  (American  Petroleum  Institute,  Recommended  Practice  for  Determining 
Permeability  of  Porous  Media,  RP  27,  Sept.  1952). 

Hydraulic  conductivity.  In  soil  mechanics  and  foundation  engineering  the 
term  permeability  is  often  applied  mistakenly  to  another  constant.  This  second 
constant  or  characteristic  typifies  not  only  the  porous  material  but  also  the  fluid 
flowing  through  it  and  includes  in  addition  to  the  permeability  as  defined  above, 
the  fluid  factors  viscosity  and  specific  weight.  In  this  report  this  product  is  termed 
"hydraulic  conductivity"  (Soil  Science  Society  of  America,  Subcommittee  on 
Permeability  and  Infiltration,  August  31,  1951).  In  terms  of  the  hydraulic  conduc¬ 
tivity  Darcy's  law  of  flow  may  be  expressed  as  follows:  The  volume-flow  per  unit 
area  is  equal  to  the  product  of  the  hydraulic  conductivity  by  the  hydraulic  gradient 
or  loss  of  head  per  unit  distance  in  the  direction  of  flow.  The  hydraulic  conduc¬ 
tivity  is  the  permeability  multiplied  by  the  specific  weight  of  fluid  and  divided  by 
its  viscosity. 

Transmissivity.  A  convenient  term  used  to  describe  the  ability  of  a  uniformly 
thick  bed  or  aquifer  to  transmit  water  is-  its  transmissivity.  This  is  the  same  as 
the  transmissibility  of  Theis,but  since  the  water  in  the  aquifer  is  transmissible  while 
the  aquifer  itself  is  transmissive,  we  prefer  the  word  "transmissivity"  to  decribe 
the  aquifer.  The  transmissivity  is  the  hydraulic  conductivity  miltiplled  by  the 
thickness  of  transmissive  formation. 

Storativity.  The  ability  of  an  aquifer  to  store  water  is  expressed  by  a 
number  called  in  this  report  its  storativity.  This  is  the  same  as  "coefficient  of 
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storage”  or  "storage  coefficient"  of  Theis.  It  is  defined  as  the  volume  of  water 
removed  from  storage  under  a  unit  surface  area  by  a  unit  decline  of  head,  or  ex¬ 
pressed  otherwise,  it  is  the  volume  of  water  stored  under  a  unit  surface  area  by 
a  unit  rise  of  head.  In  confined  aquifers  the  storativity  ranges  from  about  10”^ 
to  10“^,  in  unconfined  aquifers  from  about  .05  to  .40.  Intermediate  values  imply 
semi  -conf  i  nement . 

Pumping  level  is  the  water  level  observed  inside  a  well  when  it  is  dis¬ 
charging. 

Standing  level  is  the  water  level  observed  inside  a  well  when  its  pump 
is  idle. 

Static  level  is  the  water  level  that  would  have  been  observed  at  any 
time  in  a  well  if  it  had  never  been  pumped. 

Specific  capacity  is  the  ratio  of  the  discharge  to  the  drawdown  it  produces, 
measured  inside  the  well  and  expressed  in  gpm/ft  or  in  cfs/ft. 

Specific  drawdown  is  the  ratio  of  drawdown  to  the  discharge  that  produces 
it.  It  is  the  reciprocal  of  specific  capacity  and  may  be  expressed  in  ft/gpm  or  in 
ft/cfs. 

Residual  drawdown.  The  difference  at  any  time  between  the  standing  water 
level  observed  in  a  well  after  the  same  well  or  a  neighboring  well  has  been  shut 
drown  and  the  extrapolation  of  the  standing  water  level  from  the  next  preceding 
period  of  recovery. 

Theory  of  Well  Tests 

The  theory  of  well  tests  is  outlined  in  several  papers  that  are  listed  in  the 
Bibliography.  See  especially  Jacob[  1950]  ,  Jacob  [1947],  and  Cooper  and 

Jacob  [1946].  The  object  of  these  tests  is  three-fold - First,  to  determine  the 

formation  constants;  second,  to  determine  the  characteristics  of  the  pumping  well; 
and  third,  to  enable  predictions  of  the  future  performance  of  the  well. 

Darcy's  law.  The  law  of  horizontal  flow  of  homogeneous  fluids  through 
porous  media,  such  as  water  through  sandstone,  known  as  Darcy's  low,  is  as 
follows:  Q  =  (kA^)dp/dx,  where  Q  is  the  volume-rate  of  flow  of  water  through 
area  A  under  a  pressure  gradient  dp/dx.  The  symbolyW  stands  for  the  viscosity  of 
the  water,  and  the  symbol  k  stands  for  the  permeability  of  the  material  through 
which  the  water  is  flowing.  By  dimensional  analysis  it  may  be  seen  that  the  unit 
of  k  is  length  squared  (1^.  The  standard  unit  of  permeability  is  the  darcy. 

The  foregoing  equation  applies  only  to  horizontal  flow.  When  the  flow 
is  in  some  general  direction  the  driving  force  producing  the  flow  is  expressed  in 
terms  of  head  (h).  Also,  it  is  found  convenient  to  lump  together  the  permeability 
of  the  material  through  which  the  water  flows,  the  specific  weight  of  water  and 


also  its  viscosity^  to  give  the  following  relationship:  Q  =  KA  dh/dx,  with 
K  =  if  V/ ju  and  h  =  z  +  p/9^ 

Transmission  of  water.  Many  water-bearing  formations  in  nature  are  of 
almost  uniform  thickness,  in  which  case  it  is  found  convenient  to  lump  together 
the  factor  K,  hydraulic  conductivity,  with  the  thickness  of  the  formation  to  form 
the  coefficient  known  as  "transTiissivity" ,  as  already  stated.  Thus  Q=  Kabdh/dx 
aTdh  /dx,  where  T  ■=  Kb.  In  this  last  equation,  a  is  the  width  of  a  strip  of  the 
waterJaearing  formation  or  aquifer  at  right  angles  to  the  direction  of  the  flow  (x); 
b  is  the  average  thickness  of  the  aquifer. 

Storage  of  water.  Th e  foregoing  coefficients  relate  to  the  ability  of  an 
aquifer  to  transmit  water.  Another  important  characteristic  of  an  aquifer  is  that 
which  relates  to  its  ability  to  store  water.  This  is  called  the  "storage  coefficient" 
(Theis)  or  "storativity" ,  already  defined  above. 

Aquifers  may  be  divided  into  two  kinds--  confined  aquifers  and  unconfined 
aquifers.  vVater  may  be  stored  and  removed  from  storage  within  an  unconfined 
aquifer  by  the  rising  and  falling  of  the  water  table.  The  volume  of  water  that 
goes  into  storage  and  comes  out  of  storage  is  a  fraction  of  the  total  space  occupied 
by  the  water-bearing  material.  In  coarse-grained  materials  this  specific  volume 
may  approach  the  porosity  of  the  material.  In  fine-grained  materials,  owing  to 
capillary  action  the  volume  of  water  going  into  and  coming  out  of  storoge  may  be 
somewhat  less  than  the  actual  porosity.  The  ratio  between  the  volume  of  stored 
water  and  the  volume  of  space  it  occupies,  known  as  the  storativity  is  defined 
by  the  following  equation:  5  =  dV/Adh,  where  dV  is  the  volume  of  stored  water 
under  a  surface  area  A,  and  where  dh  is  the  differential  of  head  accompanying 
the  storage  of  a  differential  volume  dV. 

In  confined  aquifers  water  is  stored  through  the  compressibility  of  the 
water  and  the  expansibility  of  the  solid  framework  of  the  aquifer.  It  comes  out 
of  storage  through  the  expansion  of  the  water  and  the  concomitant  compression 
or  compaction  of  the  solid  framework  . 

VVel  I  Testing  -Procedures 

Interference  tests.  V/hen  a  well  that  has  been  Idle  suddenly  begins 
pumping  at  a  constant  rate  from  an  aquifer  of  uniform  thickness  and  uniform 
permeability,  the  drawdown  at  some  distance  r  away  from  that  well  and  at 
some  time  t  since  the  start  of  pumping  Is  given  approximately  by  the  following 
equation:  s  =  (2.30  Q/AnX)  log  (2.  25  Tt/Sr^,  where  S  and  T  are  the  storativity 
and  transmissivity  of  the  aquifer,  respectively,  and  Q  is  the  discharge  rate  of  the 
well.  It  is  seen  from  this  equation  that,  by  plotting  s  against  t  at  any  constant  r. 


the  points  should  fall  on  a  straight. line  on  a  semi-log  paper/  the  linear  scale  of 
which  is  used  for  plotting  drawdown  (s)  and  the  logarithmic  scale  of  which  is  used 
for  plotting  time  (t).  Moreover/  from  the  slope  of  the  straight  line  drawn  through 
such  points  plotted  on  such  a  semi-log  graph/  it  should  be  possible  to  determine 
the  transmissivity  (T)  using  this  relation:  T  =  2.  30  QZilog  t/4'vrA  s,  where  As/ A 
log  t  is  the  slope  of  the  straight  line  on  the  semi-log  graph.  For  convenience/ 

A  log  t  may  be  taken  to  be  equal  to  1  /  in  which  case  As  is  the  change  of  draw¬ 
down  over  a  ten-fold  variation  of  t/  or  over  one  "log  cycle". 

By  extending  the  straight  line  back  in  time  to  where  it  intercepts  the  line 
of  zero  drawdown  (or  time  axis)/  it  is  possible  to  determine  5  as  follows: 

5  =  2.  25  Ttg/r^,  where  t©  is  the  time  at  the  intercept  on  the  zero  drawdown 
line. 

Drawdown  tests.  The  drawdown  inside  a  pumping  well  is  made  up  of  two 
components  —  the  first/  called  "formation  loss"  /  is  the  loss  of  head  from  that  - 
great  distance  where  the  drawdown  is  negligible  up  to  the  face  of  the  well/  and 
the  second/  the  "well-loss"  is  the  loss  of  head  accompanying  the  flow  of  water 
through  the  perforations  in  the  casing  and  upward  inside  the  casing  to  the  pump 
intake.  The  formation  loss  is  given  by  an  equation  similar  to  that  given  previously 
except  that  the  effective  radius  of  the  well  (rw)  is  substituted  in  place  of  r.  Then 
writing  the  well-loss  as  CQ2  the  drawdown  in  a  pumping  well  may  be  expressed 
as  follows:  Sy^  =  (2.30  Q/4'rTT)log  (2.25  Tt/5rvy2)  +  CQ^.  This  equation  may 
be  abbreviated  as  follows:  s^  =  A(t)Q+CQ2/  where  A(t)  is  the  formation-loss 
coefficient/  which  varies  with  time/  and  where  C  is  the  well-loss  coefficient/ 
which  is  constant.  [See  Jacob/  1947.] 

It  may  be  seen  from  the  lost  two  equations  that  a  semi-log  plotting  of 
drawdown  against  time  should  again  permit  the  determination  of  T,  though  not 
of  5.  Observation  wells  at  known  distances  are  needed  to  get  S. 

Step-drawdown  tests.  By  pumping  a  well  atdifferent  rates  over  successive 
periods  of  time  and  observing  the  trend  of  drawdown  during  each  period  it  is 
possible  to  determine  the  well-loss  coefficient/  C/  and  the  formation-loss 
coefficient/  A,  which  is  a  function  of  the  duration  of  each  step.  This  kind 
of  test  is  known  as  a  "step-drawdown  test" . 

Dividing  the  foregoing  equation  for  drawdown  by  the  discharge  of  the 
well/  Syy/Q  =  A(t)  +  CQ/  the  ratio  sw/Q  being  the  specific  drawdown  of  the 
well .  A  plotting  of  specific  drawdown  against  discharge  for  several  steps  of 
a  test  should  give  astroigh'  line/  the  slope  of  which  gives  the  well-loss  coeffi¬ 
cient/  C/  and  the  intercept  of  which  gives  the  value  of  formation-loss  coefficient 
for  the  duration  of  step  that  was  used. 
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The  format  ion  “loss  coefficient  is  also  a  function  of  the  transmissivity, 
being  inversely  proportional  to  T.  Thus  it  may  be  written  A  =  A'/T,  Then, 
dividing  the  drawdown  equation  by  instead  of  by  Q,  one  may  write  s^/Q^  = 
A'/TQ  +  C.  In  this  case,  when  plotting  Sy^/Q^  against  1/QT,  the  slope  gives 
A'  and  the  intercept  gives  C. 

Not  only  may  the  foregoing  plottings  be  used  to  analyze  the  data  from 
step-drawdown  tests  of  single  wells,  but  they  may  be  used  to  analyze 
statistically  the  data  from  simple  drawdown  tests  (one-step  drawdown  tests) 
of  many  wells.  Thus  the  overage  characteristics  of  a  large  number  of  wells 
of  similar  design  may  be  determined. 

Recovery  tests.  vVhen  a  well  that  has  been  pumping  continuously  for  some 
time  is  suddenly  shut  down,  the  water  surface  inside  the  casing  begins  to  recover, 
and  continues  to  recover  for  some  time  at  a  continuously  diminishing  rate.  It 
is  just  as  though  the  well  were  to  continue  discharging  and  a  recharge  well  of 
the  same  strength  were  to  be  superimposed  on  top  of  it  at  the  time  of  shutdown. 

In  other  words,  the  position  of  the  recovering  water  surface  is  the  composite 
of  two  effects,  the  drawdown  and  the  superimposed  recovery.  This  kind  of 
test  is  known  as  a  Theis  recovery  test  (after  C.  V.  Theis  of  the  U.  i.  Geolo¬ 
gical  Survey,  see  paper  by  Theis,  Trans.  A.G.U.,  vol.  16,  p.  520,  1935). 

The  residual  drawdown  at  time  t'  after  shutdown  and  time  t  after  the  start  of 
pumping  is  given  by  this  equation:  s'  =  (2.30  QMtT T)  log  (t/t').  The  logarithms 
of  the  square  of  the  distance  and  of  the  constant  factors  would  appear  twice  in 
the  equation  for  residual  drawdown,  but  oppositely  signed,  and  therefore  they 
actually  disappear  from  the  equation. 

Tests  in  thin  unconfined  aquifers.  The  foregoing  methods  apply  to 
confined  aquifers  and  also  to  unconfined  aquifers,  provided  the  maximum  draw¬ 
downs  are  a  small  fraction  of  the  original  or  undisturbed  depth  of  flow  and 
provided  the  storativity  is  fairly  constant.  When  the  drawdown  is  a  large  fraction 
of  the  undisturbed  depth  of  flow  the  hydraulic  conductivity  is  determined  by  plotting 
values  of  h^  against  (log  r).  Knowing  K,  vajues  of  i  may  be  calculated.  As  3  is 
less  likely  to  vary  with  r  than  with  t,  calculated  values  of  5  averaged  over  the  range 
of  r  are  extrapolated  in  time  to  estimate  the  ultimate  value  of  S.  See  Jacob 
[1950,  p.  384]. 
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APPENDIX  C 

FURTHER  DETAILS  ON  WELL  TESTING  PROCEDURES 

There  follows  a  brief  summary  of  procedure  in  the  well  tests.  This  is 
given  to  supplement  the  material  found  in  Appendix  B.  That  appendix  outlines 
the  theory  of  groundwater  hydraulics  applicable  to  well  tests  and  briefly 
discusses  the  procedures  of  such  tests.  However,  it  was  felt  desirable  to 
elaborate  on  these  procedures  somewhat.  The  references  given  in  Appendix  A 
in  the  Bibliography  under  the  heading  "Well  Tests"  will  prove  useful  guides 
in  this  connection,  especially  the  text  by  Muskat  and  the  chapter  in  the  book 
"Engineering  Hydraulics"  edited  by  Hunter  Rouse.  Chapter  V  of  that  compendium, 
entitled  "Flow  of  Groundwater,"  was  written  by  the  author. 

The  drawdown-recovery  tests  should  be  run  on  as  many  wells  as  possible 
in  the  contaminated  area  that  runs  northward  from  the  northwest  boundary  of 
the  Arsenal,  especially  in  Sections  9,  IC,  15,  16,  21,  and  22  of  T2S,  R67W. 
Those  wells  in  that  area  still  equipped  with  pumps  and  prime  movers  may  be 
tested  with  little  expense.  Where  it  is  possible  to  provide  a  valve  to 
throttle  the  well,  step-drawdown  tests  may  be  run  in  addition  to  the  usual 
drawdown-recovery  test.  After  these  tests  have  been  run  and  analyzed  it 
will  be  possible  to  select  certain  ones  of  them  for  further  testing  by  the 
interference  method.  This  will  require  the  drilling  of  a  few  observation 
wells  at  each  site.  Thereafter  it  may  be  found  desirable  to  run  interfer¬ 
ence  tests  at  specially  planned  locations,  having  drilled  especially  for 
that  purpose  new  production'  wells  v/hich  might  later  be  converted  into  rec¬ 
lamation  wells  to  decontaminate  the  aquifer.  However,  this  kind  of  test 
should  be  deferred  until  all  other  possibilities  have  been  exhausted. 
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Drawdown-Recovery  Test 

This  kind  of  test  requires  no  observation  wells  but  does  require 
accurate  measurements  of  discharge  and  drawdown  in  the  pumping  well  itself, 
both  before  and  during  the  test.  See  pages  B-4  and  B-5  of  Appendix  B. 

Measurements  of  discharge.  The  discharge  should  be  measured  by  means 
of  a  circular  orifice  on  an  extension  of  the  discharge  pipe.  Where  the 
discharge  pipe  is  threaded  at  its  end  this  may  easily  be  extended.  In 
other  cases  it  may  be  necessary  to  provide  a  Dresser  coupling  followed  by 
a  valve  for  the  purpose  of  throttling  the  well  for  step-drawdown  tests,  and 
then  a  length  of  pipe  equal  to  about  10  diameters,  on  the  end  of  which  the 
orifice  is  placed.  If  it  is  not  possible  to  install  an  orifice,  the  dis¬ 
charge  perhaps  could  be  measured  by  a  sharp-edged  weir  in  a  weir  box  es¬ 
pecially  constructed  for  the  purpose  so  as  to  be  movable  from  site  to  site. 
Or  the  discharge  might  possible  be  measured  volumetr ical ly  in  a  tank  that 
could  be  transported  from  site  to  site.  However,  the  orifice  method  is  to 
be  preferred. 

Measurements  of  drawdown  inside  the  pumping  well  should  be  made  either 
by  an  electrical  probe  or  by  a  steel  tape.  In  most  instances  it  should  be 
found  possible  to  enter  the  well  through  the  bushing  in  the  pump  base  pro¬ 
vided  for  that  purpose.  A  steel  tape  that  is  somewhat  rusty  will  give  a 
better  water  line  than  a  new  tape.  The  tape  may  be  covered  with  carpenter's 
chalk  to  make  the  water  line  visible  to  the  nearest  hundreth  of  a  foot. 

Procedure.  After  a  period  of  rest,  during  which  the  standing  water 
level  is  measured  repeatedly  to  establish  its  trend,  the  discharge  of  the 
well  is  begun  and  the  trend  of  the  drawdown  during  the  period  of  pumping  is 
observed  by  frequent  measurement.  The  discharge  also  should  be  checked  fre¬ 
quently  to  see  whether  it  declines,  especially  during  the  early  part  of  the 
pumping  period. 
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After  pumping  at  a  fairly  uniform  rate  for  an  hour  or  two  the  pump 
should  be  shut  down  abruptly,  and  the  recovery  of  the  standing  water  level 
inside  the  well  should  be  measured  frequently. 

graphical,  tnninttlt*  A  semi  logarithmic  plotting  of  the  drawdown  and  of 
the  recovery  against  the  logarithm  of  the  time  elapsed  since  the  change  in 
pumpage  that  produced  that  drawdown  or  recovery  should  enable  an  estimate  of 
the  transmissivity  and  hydraulic  conductivity  of  the  aquifer.  Owing  to  the 
fact  that  the  storativity  of  the  aquifer  will  vary  in  time  and  will  vary  in 
a  different  manner  during  drawdown  than  during  recovery,  the  interpretation 
of  the  graph  may  not  be  easy. 

Step-Drawdown  Test 

Valved  discharge.  For  a  step-drawdown  test  accurate  means  of  measuring 
discharge  and  drawdown  inside  the  pumping  well  is  required,  as  also  a 
valve  for  the  purpose  of  controlling  closely  the  discharge. 

Measurements  of  drawdown  and  discharge  should  be  carried  out  the  same 
way  as  in  the  drawdown-recovery  tests. 

Procedure.  After  the  period  of  rest  during  which  the  standing  water 
level  in  the  well  is  measured  repeatedly  to  establish  its  trend,  the  dis¬ 
charge  of  the  well  is  begun  at  a  fraction  of  its  full  capacity,  say  one  third, 
and  allowed  to  continue  for  about  one  hour.  During  this  first  period  of 
pumping  the  trend  of  the  drawdown  is  observed  by  frequent  measurement.  At 
the  end  of  the  first  period  of  pumping 'the  valve  is  abruptly  opened  to  raise 
the  discharge  to  say  two  thirds  of  its  full  capacity.  Again  the  new  trend 
of  the  drawdown  is  observed  and  extrapolated  over  the  third  period,  at  the 


beginning  of  which  the  throttling  valve  would  be  opened  wide  to  allow  the 
pump  to  operated  at  full  capacity. 


Graphical  treatment.  The. data  from  the  step-drawdown  test  should  be 
plotted  on  a  graph  like  Figure  1  in  the  preliminary  report.  It  will  be 
found  that  the  formation  loss  is  curvilinear.  By  the  theory,  already 
briefly  outlined,  it  should  be  possible  to  divide  the  total  drawdown  at 
different  rates  of  pumping  into  formation  loss  and  well  loss.  However, 
because  of  the  variable  storativity  of  the  aquifer  some  special  interpre¬ 
tation  may  be  required. 


Interference  Test 

For  each  interference  test  at  least  four  observation  wells  should  be 
drilled  penetrating  sufficiently  below  the  water  table  to  allow  for  the 
drawdown  during  the  period  of  pumping.  These  may  be  located  two  each  on 
two  lines,  one  line  running  up-gradient  from  the  discharging  well  (or  down- 
gradient  if  necessary)  and  the  other  running  cross-gradient  from  the  pumping 
well..  On  each  line  there  should  be  a  well  about  50  feet  from  the  pumping 
well  and  another  about  100  feet  away.  If  it  is  desired  to  get  better  re¬ 
sults,  then  six  observation  wells  should  be  drilled  just  penetrating  the 
water  table,  three  on  each  line,  one  50  feet  away,  one  100  feet  away,  and 
one  about  200  feet  away.  In  those  places  where  the  aquifer  is  thicker,  it 
may  be  desired  to  get  still  sharper  results  by  having  pairs  of  observation 
wells  at  each  observation-well  site,  one  penetrating  just  sufficiently  far 
below  the  water  table  to  allow  for  the  drawdown  during  the  pumping  period 
and  the  other  reaching  to  the  bottom  of  the  aquifer.  However,  in  most  in¬ 
stances  this  will  not  be  necessary. 

Measurement  of  drawdown  and  discharge.  The  same  care  must  be  taken  in 

interference  tests  to  observe  accurately  the  discharge  of  the  well.  Pumping 

\ 

by  electrically  driven  turbine  is  preferred  because  of  the  constancy  of  the 
speed  and  therefore  the  near  constancy  of  the  discharge.  However,  the  dis¬ 
charge  should  be  observed  sufficiently  close  to  detect  any  declining  trend. 
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The  pumping  would  have  to  continue  24  or  48  hours  in  order  to  produce  strong 
influences  at  the  distance  out  to  the  farthest  observation  well. 

Gretphical  treatment.  The  data  from  this  test  is  best  treated  by  plot¬ 
ting  it  on  semi  logarithmic  paper,  arranging  the  horizontal  scale  or  logari¬ 
thmic  scale  to  record  time  elapsed  since  the  start  of  pumping  and  the  vertical 
scale  or  ordinate  to  record  the  drawdown  in  the  different  wells.  Theoreti¬ 
cally,  if  the  storativity  were  constant  and  the  formation  of  uniform  thickness 
and  homogeneous,  the  drawdown  curves  should  be  straight  lines  on  such  a  graph. 
However,  it  is  anticipated  that  especially  because  of  the  variableness  of  the 
storativity,  the  drawdown  data  will  plot  on  curves.  Special  care  must  be 
taken  in  their  interpretation.  Special  attention  should  be  paid  to  the  last 
paragraph  on  page  B-5  of  Appendix  B,  which  deals  with  modifications  of  the 
theory  required  for  tests  in  thin  unconfined  aquifers,  in  which  the  draw¬ 
down  is  a  sizable  fraction  of  the  undisturbed  or  initial  depth  of  flow. 

A  set  of  forms  that  have  been  found  useful  for  tabulating  the  data 
from  the  pumping  tests  should  be  used.  Space  is  provided  on  each  for  the 
calculations  leading  to  determinations  of  transmissivity  and  storativity. 
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ANALYSIS  OF  WELL  TESTS 
I  ntroduction 

During  September  and  October  of  1960  several  wells  in  Township  2  South,  Range 
67  West  were  tested.  The  locations  of  these  wells  are  given  on  Figure  D-1 .  Draw¬ 
down  recovery  tests  were  run  on  10  wells.  On  one  of  these  10  wells  a  step-drawdown 
test  was  run,  from  which  we  determine  not  only  the  hydraulic  conductivity  but  also  the 
well-loss  coefficient.  Three  of  the  ten  wells  were  used  os  discharge  wells  in  interference 
tests.  Piezometers  or  "observation  wells"  were  drilled  at  distances  of  about  50,  100  and 
200  feet  on  two  radial  lines  from  each  pumping  well.  From  these  three  tests  the  hydraulic 
conductivity  and  effective  overage  storativity  were  determined,  and  also  the  effective 
radius  of  the  pumping  well.  (For  explanation  of  symbols  used  here  see  page  D-  1  7.) 

Drawdown-Recovery  Tests 

Theory  and  procedure.  The  theory  of  drawdown  and  recovery  tests  such  as  were 
run  on  the  ten  wells  near  Rocky  Mountain  Arsenal,  are  given  in  Appendix  B,  pages 
B-4  and  B-5  and  in  the  Memorandum  of  the  Conference  of  9  August,  1960  page  4. 

Test  of  Aden  Well .  Figure  D-2  is  a  water-level  graph  for  the  drawdown-recovery 
test  of  the  Aden  Well,  Plotted  thereon  are  readings  of  depth  to  water,  in  feet,  versus 
clock  time.  Pertinent  data  on  this  test  are  os  follows: 

bate  of  test  17  September  1960. 

Started  pumping  10:37. 

Stopped  pumping  11:50, 

Measured  recovery  until  13:25. 

Discharge  declined  from  200  to  190  and  averaged  194  gpm. 

Discharge  measured  by  3-inch  orifice. 

Drilled  depth  not  known. 

Sounded  depth  24.1  ft. 

Measuring  point  1  ft  above  ground. 

Type  of  screen  or  perforated  section  not  known  , 

Lenth  of  same  not  known. 

For  other  pertinent  data  see  Table  D-1,  first  line. 

Note:  "No  water  came  up  on  first  start;  we  hod  neglected  to 
prime  pump," 

Figure  D-3  is  a  semi-log  drawdown-recovery  graph  for  the  Aden  Well,  on  which 
are  plotted  values  of  the  square  of  the  depth  of  flow  as  observed  in  the  well  against  the 


FIG.  D-2 
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logarithm  of  tho  elapsed  time.  The  square  of  the  head  difference  is  seen  to  drop 

2  2 

from  41 .9  ft  to  31 .6  ft  between  10  and  100  minutes  after  start  of  the  recovery. 

2  2 

Thus  the  change  In  h  over  one  log  cycle  is  10.3  ft  .  (See  column  [ll]  of  Table 
D-l.) 

According  to  our  approximate  theory  of  radial  flow  in  an  unconfIned  aquifer, 
developed  by  analogy  to  flow  in  a  confined  aquifer,  the  relationship  between  the 
square  of  the  depth  of  flow,  the  effective  radius  of  the  well,  the  elapsed  time,  and 
the  formation  constants  is  as  follows: 

(h  =  (2.30  Q/2  7r  K)  log  (2. 25  Kh*t/S*r3 

2 

From  the  change  in  h  over  one  log-cycle  on  the  semi-log  drawdown-recovery 
graph  the  hydraulic  conductivity  is  determined  from  the  following  relationship  which 
may  be  derived  from  the  foregoing  equation: 

K  =  2.30  Q  Adog  t)/2  rr  A(h^) 

In  this  equation  the  subsemptw  has  been  omitted  for  convenience,  as  it  has  also 
on  Figure  D-3  and  the  following  figures. 

The  measured  discharge  of  the  Aden  well  was  37,400  ft^/doy.  The  hydraulic  conductivity 
is  calculated  to  be  1,330  ft/day.  (See  columns  [10]  and  [12],  Table  D-l.) 

As  it  is  not  possible  to  determine  the  effective  average  storativity  from  a  drawdown- 
recovery  test,  and  os  it  is  not  possible  either  to  determine  the  effective  radius  of  the  well, 
these  two  quantities  are  lumped  together.  Moreover  os  the  effective  average  depth  of  flow 
is  not  known.  It  is  here  designated  h*  and  combined  with  the  other  two  unknowns.  The  logo- 
rithm  of  the  quantity  (S*  r  ^/h*  [ft])  may  be  determined  from  the  drawdown  observed  in  the 
well  at  any  particular  time,  the  hydraulic  conductivity  already  having  been  found.  For 
convenience  we  select  the  drawdown  at  100  minutes,  tabulated  in  column  [15]  of  Table  D-l , 
The  quantity  we  seek,  then,  is  given  by  the  following  formula; 

log(S*r  ^/h*[ft])  =  log  (2.25  Kt/[ft])  -  A(h^)/  A,(h^) 

2 

In  this  equation  /\(h  J  means  the  change  of  drawdown  over  one  log  cycle  on  the  semi¬ 
log  time  drawdown  graph.  In  the  case  of  the  Aden  well  this  was  10.3  feet  (Column  [1 1]). 

Refe  rring  to  Figure  D-3,  the  drawdown,  in  feet  squared  ,  after  100  minutes  of  continuous 
punfiping  was  (86.1  -  29.0)  ft^.  As  the  ratio  given  in  Column  [14]  is  5.55,  the  logarithm 
of  the  desired  quantity,  in  Column  [15],  is  2.32  -  5.55  =  -3.23  or  0.77  -  4. 

The  initial  depth  of  flow  at  the  Aden  well  is  9,3  feet,  os  given  in  Column  [16]  of 
Table  D-l. 


Tstf  gf  Dqibb^  Well.  Figure  D-4  is  a  graph  of  observed  depth  to  water  during  the  draw- 
down-recovery  test  of  the  Dobbs  well  near  Rocky  Mountain  Arsenal.  Pertinent  data  on  this 
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fest  are  as  Follows: 


Date  of  test  17  September  1960, 

Started  pumping  08:32. 

Stopped  pumping  09:46. 

Measured  recovery  until  10:15, 

Discharge  varied  from  25  to  42  and  averaged  42  gpm. 

Discharge  measured  by  3-inch  orifice. 

Drilled  depth  40  feet. 

Sounded  depth  45  Feet. 

Measuring  point  1,1  feet  above  ground. 

Type  of  screen  or  perforated  section  steel. 

Length  of  same  4  inches  to  6  inches. 

For  other  pertinent  data  see  Table  D-1 ,  second  line. 

Notes:  "Well  not  used  for  one  year.  Owner  claimed  discharge 
about  400  gpm.  Changed  to  small  orifice  at  08:34," 

Figure  D-5  is  a  semi-log  drawdown-recovery  graph  of  the  Dobbs  well.  The  points 
plotted  there  os  circles  represent  the  drawdown  and  those  plotted  as  dotS/  the  recovery.  The 
recovery  was  measured  up  from  the  extrapolated  drawdown  and  referred  to  a  new  beginning  in 
time,  namely  the  time  of  shutdown.  In  an  ideal  situation  this  transformation  should  bring  the 
two  sets  of  data  into  sup*erpositiois.  The  spread  of  the  data  is  attributable  to  inhomogeneities 
in  the  formation  and  also  to  the  variability  of  the  storativity. 

In  this  case  the  change  of  h^  over  one  log-cycle  is  only  about  2,5  ft  and  the  hydraulic 
conductivity  about  1,190  feet  per  day.  The  log  of  ^S*r^  /h*[ft])  is  0,94  -  6,  (See  Columns 
[12]  to  [15]  of  Table  D-1.) 

Ta«t  of  Marty  WelL  Figure  D-6  is  a  drawdown-recovery  graph  for  the  test  of  the  Marty 
well  near  Rocky  Mountain  Arsenal.  Observations  of  depths  to  water,  in  feet,  are  plotted  there 
against  the  clock  time.  Pertinent  data  on  this  test  are  os  follows: 

Date  of  test  21  September  1960.  _ 

Started  pumping  08:55. 

Stopped  pumping  10:00. 

Measured  recovery  until  15;  16. 

Discharge  declined  from  886  to  732  and  averaged  765  gpm. 

Discharge  measured  by  5-inch  orifice. 

Drilled  depth  52  feet. 

SourKled  depth  not  given. 

Measuring  point  Flush  with  grourKi  surface. 

Type  of  screen  or  perforated  section  concrete. 
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Length  of  well  screen  not  known. 

For  other  pertinent  data  see  Table  D-1,  third  line. 

Figure  E-7  is  a  semi-log  drawdown-recovery  graph  for  the  Marty  well.  The  circles 

2 

represent  calculated  values  of  h  determined  from  observations  of  depth  to  water  during 

the  pumping  and  subsequent  recovery.  The  dots  represent  the  recovery  replotted  on  the 

logarithmic  time  scale,  taking  the  time  of  shutdown  as  the  new  origin. 

2 

The  change  of  h  over  one  log-cycle  is  18  feet, -and  the  hydraulic  conductivity  is 
calculated  to  be  about  3,000  ft/day.  The  "drawdown"  at  100  minutes  is  377  -  76  or  301  ft^. 
Accordingly  the  log  of  (S*r  /h  [ft])  =  .97  -  15.  The  reason  for  the  high  negative  value 
here  is  that  the  well  is  strong  and  the  drawdown  is  large  in  comparison  to  the  initial  depth  of 
flow.  This  entails  high  convergence  losses  as  the  water  decends  toward  the  well.  The  initial 
depth  of  flow  in  this  cose  was  19.4  feet. 

Test  of  Masunogg  Well .  Figure  D-8  is  a  drawdown-recovery  graph  of  the  test  of  the 
Masunaga  well  near  Rocky  Mountain  Arsenal.  The  well  was  first  pumped  for  about  4  minutes 
starting  at  12:05,  when  it  was  found  necessary  to  reduce  the  size  of  the  orifice.  After  wait¬ 
ing  several  minutes,  pumping  was  resumed.  It  is  believed  that  this  short  initial  pumping 
period  after  the  false  start  hod  little  effect  on  the  subsequent  test  and  accordingly  corrections 
were  not  made  for  it.  Pertinent  data  on  this  test  are  os  follows: 

Date  of  test  20  September  1960. 

Started  pumping  12:41. 

Stopped  pumping  13:40. 

Measured  recovery  until  14:15.  (Two  measurements  were 
made  the  following  day.) 

Discharge  declined  from  183  to  169  and  averaged  176  gpm. 

Discharge  measured  by  3-inch  orifice. 

Drilled  depth  39  feet. 

Sounded  depth  not  given. 

Measuring  point  0.9  foot  above  ground. 

Type  of  screen  or  perforated  section,  concrete. 

Length  of  same,  12  feet  (?). 

For  other  pertinent  data  see  Table  D-1,  fourth  line. 

Notes:  "In  this  test  we  observed  the  well  pumped  and  a  well 

which  was  [some  distance]  away.  The  well  reading  in 
the  remarics  column  [of  the  data  sheet]  has  a  correction 
of  +0,25  ft  to  each  reading.  Well  1355  [Masunaga  Well] 
is  connected  to  two  wells  by  siphons.  When  we  pumped, 

we  took  readings  on  Well  1355  and  Well  A.  I  was  not 
aware  of  Well  B  until  pumping  hod  been  completed." 


FIG.  D-8 
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Figure  D-9  is  a  semi-log  drawdown-recovery  graph  for  the  Masunaga  Well  ,  Again 

2 

the  circles  represent  calculated  values  of  h  plotted  against  the  logarithm  of  time  during 
the  drawdown  and  subsequent  recovery.  The  dots  represent  the  recovery  referred  back  to 
the  new  time  origin  and  expressed  as  a  differential  of  h  measured  up  from  the  extrapolated 
drawdown  curve. 

The  change  of  h'^  over  one  log  cycle  is  in  this  case  12  feet,  and  the  hydraulic  conduc¬ 
tivity  is  calculated  to  be  1,040  ft/day.  The  drawdown  at  100  minutes  was  102,8  -  39.8  or 
about  63  ft^.  The  log  of  (S*r^^/h*[ft])  is  0,96  -  4,  In  this  case  the  initial  depth  of  flow  was 
only  10,1  ft. 

Test  of  Motsumoto  Well.  Figure  D-1G  is  a  drawdown-recovery  graph  for  the  test  of  the 
Matsumoto  well  near  Rocky  Mountain  Arsenal.  Pertinent  data  on  this  test  are  as  follows: 


Date  of  test  22  September  1960. 

Started  pumping  08:00  , 

Stopped  pumping  09:05. 

Mveasured  recovery  until  13:35,  (One  measurement  following  morning.) 
Discharge  declined  from  357  to  326  and  averaged  336  gpm. 

Discharge  measured  by  4-inch  orifice. 

Drilled  depth  not  given. 

Sounded  depth  36,9  feet. 

Measuring  point  2.75  feet  above  ground. 

Type  of  screen  or  perforated  section,  concrete. 

Length  of  same,  6  feet. 

For  other  pertinent  data  see  Table  D-1 ,  fifth  line. 


Figure  D-1 1  is  a  semi-log  drawdown-recovery  graph  for  the  Matsumoto  veil .  The 
change  of  h^  over  one  log— cycle  is  in  this  cose  9.0  ft,  and  the  hydraulic  conductivity  is 
calculated  to  be  2,640  ft/day.  The  drawdown  at  100  minutes,  expressed  in  terms  of  h  is 
131.8  -  44.0,  or  87.8  ft^.  The  log  of  (S*r^^/h*[ft])  is  0.85  -  8.  In  this  case  the  initial 
depth  of  flow  was  11.5  feet. 

Test  of  Matsumoto  Well  2.  Figure  D-12  is  d  drawdown-recovery  graph  for  the  test 
of  Matsumoto  Well  2  near  Rocky  Mountain  Arsenal.  Pertinent  data  for  this  test  are  os  follows: 


Date  of  test  22  September  1960. 

Started  pumping  11:55. 

Stopped  pumping  13:00. 

Measured  recovery  until  13:50  (one  measurement  following  morning). 
Discharge  declined  from  447  to  368  and  averaged  384  gpm. 
Discharge  measured  by  4- inch  orifice. 

Drilled  depth  45  feet. 


22  September  1960 
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Sounded  depth  36.9  feet. 

Measuring  point  1 ,0  feet  above  ground. 

Type  of  screen  or  perforated  section,  concrete. 

Lenth  of  some,  6  feet. 

For  other  pertinent  data  see  Table  D-1 ,  sixth  line. 

Figure  D-1 3  is  a  semi-log  drawdown-recovery  graph  for  Matsumoto  Well  2.  The 
same  convention  is  used  as  before  in  the  plotted  points.  The  change  of  h  over  one  log- 
cycle  is  10  ft,  and  the  hydraulic  conductivity  is  calculated  to  be  2,720  ft/doy.  The  draw¬ 
down  at  100  minutes  expressed  in  terms  of  h^  is  201  i9  -  54.0,  or  147.9  ft  .  The  log  of 
(S*r^^/h*[ft] )  is  0.83  -  13.  The  reason  for  the  high  negative  logarithm  here  is  that  the 
well  is  strong  and  the  drawdown  a  larger  fraction  of  the  initial  depth  of  flow  than  ordinarily. 
In  this  cose  the  initial  depth  of  flow  was  14.2  feet. 

Test  of  Monson  Wei  I .  Figure  D-14  is  a  drawdown-recovery  graph  for  the  test  of  the 
Monson  well  near  Rocky  Mountain  Arserwl .  Pertinent  data  for  this  test  are  os  follows: 

Date  of  test  20  September  1960. 

Started  pumping  08:25. 

Stopped  pumping  09:30. 

Measured  recovery  until  11:41  (one  measurement  afternoon  of  following  day). 

Discharge  declined  from  400  to  377  and  averaged  386  ^m. 

Discharge  measured  by  4-inch  orifice. 

Drilled  depth  40  feet,  plus  or  minus. 

Sounded  depth  not  given. 

Measuring  point  1 .7  feet  above  ground. 

Type  of  screen  not  given. 

Lenth  of  same,  not  given. 

For  other  pertinent  data  see  Table  D-1,  seventh  line. 

Notes:  "Well  in  operation  morning  of  19  September  1960.  Owners' 
well  logs  have  been  sent  to  State  Engineer's  office." 

Figure  D-15  is  a  semi-log  drawdown-recovery  graph  for  the  Monson  well .  The  change 
of  h^  over  one  log-cycle  is  12.5  feet,  and  the  hydraulic  conductivity  is  calculated  to  be 
2,180  ft/doy.  The  drawdown  at  100  minutes  is  148.9  -  76.0,  or  about  72,9  ft^.  The  log  of 
(S*r^^/h*[ft])  is  0,70  -  4,  In  this  case  the  initial  depth  of  flow  was  12,2  feet. 

Test  of  Myers  Well .  Figure  D-16  is  a  drawdown-recovery  graph  for  the  test  of  the 
Myers  well  near  Rocky  Mountain  Arsenal .  Pertinent  data  on  this  test  are  as  follows: 

Date  of  test  21  September  1960. 

Started  pumping  13:15. 


FIG.  D-14 


20  September  1960 


FIG.  D-15 
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Stopped  pumping  14:10. 

Measured  recovery  until  14:50. 

Discharge  declined  from  54  to  42  and  averaged  46  gpm. 

Discharge  measured  by  3-inch  orifice. 

Drilled  depth  not  given. 

Sounded  depth  39.5  feet. 

Measuring  point  flush  with  ground  surface. 

Type  of  screen  not  given. 

Length  of  screen^  4  feet  (?). 

For  other  pertinent  data  see  Table  D-1 ,  eighth  line. 

Notes:  "The  discharge  figures  are  estimated.  I  plotted  head  against 
drawdown  and  extrapolated  to  zero." 

Figure  D-17  is  a  semi-log  drawdown-recovery  graph  for  the  Myers  well .  The  change 
of  h^  over  one  log-cycle  is  less  than  2.5  ft,  and  the  hydraulic  conductivity  is  calculated 
to  be  more  than  1/300  ft/dw.  The  drawdown  at  100  minutes  is  12.25  +  .05,  or  about 
12.3  ft^.  The  log  of  (S*r^^/h*[ft])  is  0.41  -  3.  In  this  case  the  initial  depth  of  flow  was 
only  3.5  ft. 

Test  of  Wolpert  Well.  Figure  D-18  is  a  drawdown-recovery  graph  for  the  test  of 
the  Wolpert  well  near  Rocky  Mountain  Arsenal.  Pertinent  data  for  this  test  are  as  follows: 

Date  of  test  23  September  1960. 

Started  pumping  08:25. 

Stopped  pumping  09:31 . 

Measured  recovery  until  10:30. 

Discharge  declined  from  282  to  274  and  averaged  275  gpm.  (Initial 
discharge  262  gsm.) 

Discharge  measured  by  5-inch  and  3-lnch  orifices. 

Drilled  depth  not  given. 

Sounded  depth  18.1  feet. 

Measuring  point  flush  with  grouncf  surface. 

Type  of  screen  of  perforated  section/  concrete. 

Lenth  of  same/  not  given. 

For  other  pertinent  data  see  Table  D-1/  tenth  line. 

Note:  “Loss  1  g>m." 

Summary  of  Tests.  Table  D-1  gives  a  summary  of  the  ten  drawdown-recovery  tests  run 
on  wells  near  Rocky  Mountain  Arsenal.  Included  among  these  ten  is  the  Powers  well/  which 
was  given  a  step-drawdown  test.  It  is  reported  on  more  fully  in  the  next  section  of  this 
Appendix. 
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Column  [l]  of  Table  D-1  gives  the  name  of  the  owner  of  the  well  and  column  [2] 
the  location  by  section  number  and  fraction  thereof  within  T25  R67W.  Column  [3]  gives 
the  depth  of  the  well^  in  feet.  Those  depths  shown  to  the  nearest  foot  are  drilled  depths, 
while  those  shown  to  tenths  or  hundreths  are  depths  that  were  sounded  during  the  tests. 

Columns  [4]  and  [5]  give  the  diameter  of  the  casing  in  inches  and  the  kind  of  material 
of  which  constructed,  respectively.  Columrw  [6]  and  [7]  give  the  diameter  of  the  pump  in 
inches  and  its  power  in  horsepower,  respectively. 

In  inventorying  the  wells  each  owner  was  asked  the  estimated  discharge  of  his  well. 

In  every  cose  but  one  the  owners  overestimated  the  discharge,  sometimes  by  two  or  three 
times.  (Compare  Columns  [8]  and  [9]  of  Table  D-2.) 

The  discharge  measured  by  pipe  orifice  during  the  test  is  expressed  in  gpm  in  Column 
[9]  and  in  thousarxls  of  ft^/day  in  Column  [10].  The  figures  in  Columns  [l  1]  through  [16] 
hove  already  been  discussed  in  the  foregoing  discussion  of  the  individual  tests. 

The  overage  depth  of  the  ten  wells  tested  is  37.2  feet.  The  most  common  casing 
diameter  48  inches  arx!  the  most  common  material  concrete.  The  average  initial  depth  to 
flow  is  11.6  feet.  The  average  discharge  287  gpm,  or  55,200  ft  /day.  These  average 
figures  are  given  at  the  bottom  of  the  table.  Also  given  there  are  averages  calculated 
after  excluding  the  two  starred  numbers  in  each  column,  namely  those  for  the  Marty  well  and 
the  Matsumoto  No.  2  well.  These  two  wells  are  excluded  because  of  their  inordinately  high 
discharge  and  strong  drawdown  in  comparison  to  the  initial  depth  of  flow.  The  remaining 
eight  wells  have  an  overage  discharge  of  215  gpm  or  41,400  ft^/doy.  For  these  eight  tests 
the  overage  hydraulic  conductivity  was  about  1,500  ft/doy.  The  average  log  of  (SV^Vh^ft]) 
for  the  eight  wells  is  0.72  -  5  and  the  average  initial  depth  of  flow  10.3  feet. 

In  order  to  calculate  the  future  drawdown  around  similar  wells  it  may  not  be  necessary 
to  separate  out  the  effective  average  storatlvity,  the  effective  average  depth  to  flow,  ard 
the  effective  radius  of  the  well.  Because  of  the  crudity  of  the  theory  there  may  be  uncertain¬ 
ties  both  in  S*  and  in  h*.  Variations  in  one  or  both  may  be  reflected  in  unreal  values  of  r^ 
if  attempts  are  made  to  calculate  that  magnitude.  For  our  purpose,  we  may  use  0.72  -  5  or 
-4.28  as  the  logarithm  of  this  combined  quantity  in  predicting  the  behavior  of  similar  wells 
over  long  continued  pumping. 
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Step-Drawdown  Test  of  Powers  Well 


Theory  and  procedure  „  The  theory  and  procedure  of  step-drawdown  tests  is  given 

on  page  B-4  of  Appendix  B  arxl  on  page  5  of  Memorandum  of  Conference  of  9  August  1960, 

Pertinent  Dota  on  test.  Figure  D-20  is  a  graph  of  the  step-drawdown  test  of  the 

Powers  well  near  Rocky  Mountain  Arsenal .  At  the  bottom  of  the  graph  calculated  values 

2  2 

of  h  ^  in  ft  ,  are  plotted  against  the  clock  time.  At  the  top  of  the  graph  discharge  reod- 
ingS/  in  gpm,  are  plotted  on  the  same  time  scale.  There  are  six  steps  altogether  in  the  test. 

The  well  was  started  at  full  capacity  and  pumped  for  about  20  minutes.  After  being  allowed  to 
recover  for  about  45  minutes  it  was  started  up  again  at  a  fraction  of  full  capacity,  os  indicated 
on  the  graph.  Finally,  after  increasing  the  discharge  in  two  subsequent  steps,  the  well  was 
shut  off  and  allowed  to  fully  recover.  Pertinent  data  on  this  test  are  os  follows: 


Date  of  test  19  September  1960. 

Started  pumping  12:01. 

Stopped  pumping  (final)  15:20. 

Measured  recovery  until  16:38  (one  measurement  following  morning). 

Discharge  varied  from  step  to  step. 

Discharge  measured  by  5-inch  orifice. 

Drilled  depth  40  feet. 

Sounded  depth  not  given. 

Measuring  point  1 . 1  feet  above  ground. 

Type  of  screen  or  perforated  section,  concrete  with  l/4-4nch  slots,  V-shaped. 

Length  of  same,  10  feet. 

For  other  pertinent  data  see  Table  D-1,  ninth  line;  also  Table  D-2. 

Notes;  "Well  not  used  for  four  years.  The  owner  had  estimated  800  gpm  +. 
We  installed  a  5-inch  orifice.  Pumped  for  20  minutes,  stopped 
pump,  and  allowed  well  to  recover.  Changed  to  3-inch  orifice  and 
started  pumping  at  13:05  for  a  step-drawdown  [test].  Had  to  stop 
pumping  at  15:20,  beginning  to  flood  crops." 


Figure  D-21  is  a  semi-log  incremental  drawdown  graph  for  the  step-drawdown  test  of 
the  Powers  well .  Plotted  thereon  are  calculated  values  of  the  increment  of  h  divided  by  the 
corresponding  increment  of  discharge,  but  multiplied  by  the  increment  of  discharge  during  the 
first  step  in  order  to  keep  the  dimensions  [ft  ].  The  data  for  the  different  steps  are  distinguished 
by  distinctive  symbols  and  the  number  of  the  step  is  indicated  on  the  first  point  for  each  line  of 
data.  It  is  seen  that  the  curves  for  the  last  four  steps  agree  fairly  well .  According  to  the  theory 
of  step  drawdown  tests  in  confined  systems,  which  we  have  modified  to  apply  to  the  unconfined 
or  nonlinear  systems,  if  there  were  no  well  loss,  these  graphs  should  coincide  ,  except  for 


SEMI -LOG  INCREMENTAL  DRAW: 

FOR  POWERS  WELL 
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disturbances  caused  by  inhomogenities  and  variations  in  permeability  or  storativity. 

The  equation  for  the  drawdown  during  any  one  step  of  the  test  is  as  follows; 

ACh  ^)/h*  =  (2.30  AQ/2trKh*)log(2.25  Kh*t/Sv3  +  ^'Q^/h  ^ 

w  w 

In  this  cose  we  have  divided  the  square  of  the  head  or  depth  of  flow  measured  in  the 
well  by  the  effective  average  depth  of  flow  in  order  to  express  drawdown  in  terms  of  length 
instead  of  length  squared  os  before.  This  is  necessary  in  order  to  permit  addition  of  the 
second  term  in  the  foregoing  equation^  which  represents  the  well  loss,  the  first  term  then 
representing  the  formation  loss.  The  well  loss  is  seen  to  be  proportional  to  some  coefficient 
C  and  to  the  square  of  the  discharge.  Also  it  is  assumed  to  be  inversely  proportional  to  the 
square  of  the  head  measured  in  the  well  (h3. 

Dividing  the  differential  of  h  by  the  increment  of  discharge  that  produced  it  and 
multiplying  by  the  increment  of  discharge  during  the  first  period  we  get  the  following  relation: 

A(h3  AQ/  AQ  =  (2.30  AQ/2rr  K)log(2.25  Kh*t/S*r3 

Now  this  last  equation  is  written  assuming  C  =  0,  or  in  other  words  assuming  that  the 
vMcll  loss  is  negligable.  (In  the  case  of  the  Powers  well  we  found  this  to  be  so,  and  this  would 
seem  to  be  reasonable  in  view  of  the  fact  that  the  maximum  discharge  of  the  well  was  only 
about  293  gpm.) 

Failure  of  the  data  for  the  first  and  second  steps  of  the  drawdown  curve  to  conform  to 
the  data  for  the  other  steps  on  Figure  D-21  is  probably  attributable  to  the  fact  that  we  hove 
underestimated  the  discharge  during  the  first  period  of  the  test.  As  indicated  in  the  notes  under 
the  foregoing  tabulation  of  data  on  this  test,  it  was  found  necessary  to  interrupt  the  pumping  to 
change  from  o  5-  inch  to  a  3-  inch  orifice  because  the  owner  had  overctsti mated  the  capacity 
of  his  well.  It  is  not  unlikely  that  the  discharge  during  the  first  period  of  the  test  was  even 
greater  than  plotted  on  the  Figure  D-20.  Actually  in  (Staining  the  bottom  two  curves  on 
Figure  D-21  we  used  a  discharge  increment  of  268  gpm,  and  the  same  os  for  steps  5  arxi  6. 

To  bring. the  data  into  agreement  it  would  be  necessary  to  assume  that  the  discharge  increment 
during  the  first  and  second  steps  was  about  28  percent  greater  than  during  the  fifth  and  sixth 
steps,  or  in  orther  words  about  340  gpm  instead  of  280  to  293  gpm  as  recorded  using  the  5-inch 
orifice,  which  proved  |o  be  too  large. 

The  change  of  h  over  one  log- cycle  is  in  this  case  18  feet,  arxl  the  hydraulic  conduc¬ 
tivity  is  calculated  to  be  1 ,040  ft/day.  The  drawdown  after  100  minutes  is  about  1 18  feet 
and  the  log  of  ^  about  0.66  -  5.  The  initial  depth  to  flow  was  1 1 .8  feet. 

Table  D-2  summarizes  the  data  from  the  step-drawdown  test  of  the  Powers  well.  It 

gives  increments  of  drawdown  and  discharge  measured  20  minutes  after  each  change  in  the 

rate  of  discharge.  The  first  column  gives  the  number  of  the  step  and  the  second  and  third 

2 

colurrms  give  the  increments  of  discharge,  in  gpm  and  ft  /day  respectively.  The  fourth  and 


TABLE  D-2 


DATA  FROM  STEP-DRAWDOWN  TEST  OF  POWERS  WELL 
(increments  of  drawdown  measured  at  20  minutes.) 


[1] 

[2] 

[3] 

[4] 

[5] 

[6] 

[7] 

Increment  of 

Discharge 

_ lal _ 

Increment  of 
drawdown  ^ 

lAhVfft®] 

Sum  of  drawdown 
increments, 
ft^ 

step 

gpm 

f  t  V day 

gpm 

ftVday 

1 

+268 

+51 ,600 

268 

51 ,600 

121.2 

121.2 

2 

-268 

-51 ,600 

268 

51 ,600 

121.6 

121.6 

3 

+195 

+37,500 

195 

37,500 

70.5 

70.5 

4 

+  22 

+  4,300 

217 

41 ,800 

7.8 

78.3 

5 

+  51 

+  9,800 

268 

51 ,600 

20.4 

98.7 

6 

-268 

-51 ,600 

268 

51 ,600 

93.6 

93.6 

D-11 
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Fifth  columns  give  the  actual  discharge  or  sum  of  increments  of  discharge,  in  gpm  and  ft  /day 

respectively.  The  sixth  column  of  Table  D-2  gives  the  increment  of  drawdown  For  each  step 

2  2 

of  the  test,  expressed  as  Ah^  in  ft^.  Lcstly,  the  seventh  column  of  the  table  gives  the  sum 

2 

of  the  drawdown  increments  in  Ft  . 

Figure  D-22  Is  a  specific-capacity  diagram  For  the  Powers  well,  on  which  are  plotted 
for  each  of  the  six  steps  of  the  test  the  sum  of  the  drawdown  increments  in  square  Feet  against 
the  sum  of  the  discharge  increments  in  thousands  of  cubic  feet  per  day.  it  is  seen  that  the 
points  For  the  lost  Four  steps  of  the  test  Fall  Fairly  close  to  a  straight  line  passing  through 
the  origin,  whereas  the  points  for  the  First  and  second  steps  of  the  test  lie  somewhat  below 
that  line.  The  reasons  for  this  are  os  already  stated,  namely  that  the  discharge  during  the 
first  period  of  the  test  undoubtedly  hod  been  underestimated. 

We  may  conclude  from  this  graph  or  specific-capacity  diagram  that.  For  practical 
purposes,  over  the  range  of  discharge  involved  here,  when  dealing  with  wells  similar  in 
design  to  the  Powers  well,  the  well  loss  may  be  neglected.  Departures  From  linearity  of 
flow  in  the  vicinity  of  the  well  and  convergence  losses  as  the  water  descerxisto  the  well  may 
be  taken  care  of  by  high  negative  values  of  the  log  of  the  effective  radius  of  the  well. 


FIG.  D-22 
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Interference  Test  of  Powers  Well 


Theory  and  procedure.  The  theory  and  procedure  of  Interference  tests  of  the  kind 
run  on  the  Powers  well  is  given  on  page  B-3  of  Appendix  B  and  on  page  6  of  the  Memorandum 
of  Conference  of  9  August  i960.  Figure  D-23a  isa  map  showing  the  lay-out  for  the  interfer¬ 
ence  test  of  the  Powers  well.  Indicated  thereon  are  the  locations,  well-head  elevations  and 
depths  for  each  of  the  four  piezometers,  for  Well  1072,  Well  A  and  the  Powers  Well  itself. 

Conduct  of  test.  Well  A  is  connected  to  the  Powers  well  through  a  buried  siphon.  The 
influence  of  the  pumping  from  the  Powers  well  on  the  water  level  in  Well  A  may  be  seen  on 
Figure  D-28.  No  attempt  was  made  to  correct  for  this  effect,  however  its  neglect  is  not  thought 
to  be  of  serious  consequence. 

Before  the  well  was  pumped  the  trend  of  water  levels  was  established  in  all  the  neighbor¬ 
ing  wells  by  daily  measurements  from  5  October  through  8  October.  The  pump  was  started 
at  08:23  on  9  October  1960.  Pumping  continued  at  the  overage  rate  of  264  gpm  or  about 

3 

50,900  ft  /day  until  08:33  on  12  October  1960.  Observations  of  recovery  were  continued 
with  regularity  until  about  midnight  that  night.  Thus  the  length  of  the  pumping  period  was 
about  4,300  minutes,  and  of  the  recovery  period  about  950  minutes.  Further  infrequent  checks 
on  the  recovery  were  made  on  the  succeeding  days,  however. 

The  drawdown  and  recovery  observed  in  the  Powers  well  and  in  surrounding  wells  during 
the  test  from  October  9  to  13  are  plotted  on  Figures  D-24  through  D-30,  aS  Follows: 

D-24  Drawdown  and  Recovery  Graph  of  Powers  Well 
D-25  Graph  Showing  Interference  of  Powers  Well  on  PI 
D-26  P2 

D-27  P3 

D-28  Well  A 

D-29  P4 

D-30  Well  1072 

Analysis  of  data.  Shown  on  each  of  these  graphs  but  two,  as  a  dashed  inclined  line 

near  the  top  of  the  graph,  is  the  extrapolated  initial  trend  of  the  water  level.  The  drawdown 

is  obtained  by  taking  differences  at  anytime  between  the  values  of  h  from  tKs  line  and  the 

2 

values  of  h  read  from  the  curve  drawn  through  the  overage  of  the  plotted  points.  Similarly 

2 

the  recovery  is  taken  as  the  difference  at  any  time  between  the  h  read  from  the  recovery 
curve  traced  through  the  plotted  points  and  the  extension  of  the  drawdown  curve  at  the  same 
time.  These  differences  in  h  ,  in  square  ft,  are  plotted  on  the  semi-log  graphs  in  Figures 
D-31  through  D-35.  These  graphs  are  os  follows: 

D-31  Semi-log  Drawdown  Graph  of  Powers  Well 
D-32  Drawdown 

D-33  Recovery 


PI,  P2  and  P3 
PI .  P2  and  P3 


FIG.  D-23a 


MAP  OF  PUMPING  TEST  LAYOUT  FOR  POWERS  WELL 
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FIG.  D-32 


fig.  D-33 


fig.  D-34 
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D-34  Semi-log  Drawdown  Graph  of  P4/  Well  1072 
D-35  Recovery  P4 


On  Figure  D-31  the  differences  in  h  are  plotted  as  ordinates  against  the  time  in 
minutes  on  the  logarithmic  scaleo  In  all  the  other  graphs  the  logarithmic  scale  express 
the  ratio  (t/r  )  in  [min/ft  ],  This  composite  coordinate  is  used  to  bring  the  data  from 
the  different  wells  on  a  given  line  into  superposition,  or  os  nearly  so  as  the  inhomogenei¬ 
ties  and  variations  in  storativity  permit. 

The  points  on  the  semi-log  drawdown  and  recovery  graphs  plot  very  closely  to 
straight  lines  after  the  initial  period  of  expansion  of  the  transient.  Straight  lines  are 
drawn  through  the  data,  in  some  cases  averaging  the  data  for  two  observation  wells.  From 
the  slopes  of  these  straight  lines  it  is  possible  to  estimate  the  hydraulic  conductivity,  and 
from  their  intercept  on  the  zero  drawdown  line  the  ratio  of  effective  average  depth  of 
flow  to  effective  overage  storativity  may  be  determined. 

The  equation  for  the  drawdown  expressed  as  a  change  of  h  is  as  follows: 

A(h^)  =  (2.30  Q/2  tt  K)log(2.25  Kh*t/S*r^) 


The  hydraulic  conductivity  is  calculated  from  the  following  relation: 

K  =  2.30  Q  Adog  t)/2n  A(h^) 

At  the  intercept,  log  =  0,  and  2.25(Kh*/S*)  (t/r^)^=  1. 

Solving  this  for  the  desired  ratio  of  effective  overage  depth  of  flow  to  effective 
average  storativity,  one  gets 

hVS*  =  l/[2.25  K(t/r\] 

Table  D-3  gives  a  summary  of  the  data  from  the  interference  test  of  the  Powers 
well.  The  first  column  gives  the  name  of  the  observation  well,  including  the  Powers 
well  itself  as  one  of  these.  The  second  column  gives  the  distance  in  feet  from  each  ob¬ 
servation  well  or  piezometer  to  the  pumping  well.  No  radius  is  given  for  the  Powers 
well  itself  as  its  effective  radius  is  os  yet  undetermined. 

Column  [3]  gives  the  change  over  one  log-cycle  of  the  absolute  value  of  differen- 
ces  in  h'^,  in  fr  .  Drawdown  is  distinguished  from  recovery  by  "d",  os  against  "r*'. 

Values  of  calculated  hydraulic  conductivity  are  given  in  Column  [4]  in  ft/day. 
These  range  from  1,020  to  1,170  and  overage  1,090, 

The  intercepts  on  the  zero-drawdown  line  are  given  in  column  [5],  and  from  these 
values  of  the  ratio  (h*/S*)  given  in  Column  [6]  were  obtained.  Those  from  the  drawdown 
are  listed  on  the  left  and  those  from  the  recovery  on  the  right.  The  range  here  is  large, 
but  it  must  be  remembered  that  it  is  the  logarithm  of  this  ratio  that  is  involved  in  our 
calculations  or  predictions.  The  average  of  all  six  calculated  ratios  is  1 16  ft.  As  the 


TABLE  D-3 


DATA  FROM  INTERFERENCE  TEST  OF  POWERS  WELL 


[1] 

[2] 

[3] 

.  [4] 

[5] 

[6] 

Observation 

well 

Distance 

r/[ft] 

Chamge  over^* 
one  log-cycle  of 

lA(h")l/[ft^] 

Hydraulic 

conductivity 

K/[ft/day] 

Intercept 

h*/S*^^ 

ifti 

Powers 

— 

16.5 

1,130 

— 

PI,  2 

52.3,  99.6 

18.0  d 

1,040 

0.0042 

147 

16.0  r 

1,170 

.0025 

219 

P3 

195.3 

18.0  d 

1,040 

.0053 

116 

17.8  r 

1,050 

.039 

16 

P4 

131 .9 

18.3  d 

1,020 

.0081 

77 

16.1  r 

1,160 

.0045 

123 

Average 

1,090 

113d  119r 

a.) 

The  d  means  drawdown;  the  r  means  recovery. 
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effective  average  depth  to  flow  cannot  exceed  the  initial  depth  to  flow,  namely  11,8  ft,  we 
will  certainly  overestimate  the  overage  storativity  for  the  period  of  pumping  if  we  use  11,8  ft 
for  h*.  Doing  this,  we  get  the  following  result  S*  =  11.8  ft/1 16  ft  =  0. 102. 

Interpretation,  It  would  appear  that  the  overage  storativity  for  the  3-day  interference 
test  is  about  10  percent.  The  storativity  remained  essentially  constant  throughout  this  period, 
indicating  by  its  constancy  and  low  value  that  there  may  be  fine-grain  sediments  in  the  upper 
part  of  the  aquifer  overlying  the  more  permeable  gravel  that  cause  capillary  retention.  It 
is  not  known  how  extensive  this  condition  is.  If,  as  in  this  case,  up  to  as  much  as  20  percent 
of  the  total  volume  may  be  occupied  by  capillary  water,  there  may  be  contaminants  or 
toxicants  held  above  the  water  table  that  hove  not  been  inventoried  by  our  volumetric  cal¬ 
culations. 
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Interference  Test  of  Aden  vVell 


Theory  and  procedure.  Some  os  before. 

Conduct  of  the  test.  The  Aden  well  was  pumped  continuously  for  about  4,300  minutes, 
beginning  16  November  1960  and  ending  19  November  1960.  The  recovery  was  measured  at 
frequent  intervals  for  another  1800  minutes  and  daily  thereafter  for  several  days. 

Three  piezometers  were  drilled  on  each  of  two  lines,  one  on  an  azimuth  of  44°  29'  and 
the  other  on  an  azimuth  of  135°  10'.  Their  respective'distances  from  the  Aden  well  are  os 
follows: 


Ml! 


me 


SE  line 


PI 

50.2 


P2 

100.2 


P3 

199.3 


P4 

51.3 


P5 

98.9 


P6 

200.6  ft 


Water  level  readings  made  in  the  Aden  well  and  in  the  six  piezometers,  converted 

2  2 

to  read  in  terms  of  the  square  of  the  depth  of  flow  (h  ),  in  ft  ,  are  plotted  on  Figures  D-36a 
through  D-38b. 

Except  for  the  35-minutes  failure  in  the  pump  operation  at  about  2150  minutes  after  the 
start  of  pumping,  the  pumping  was  continous. 

Analysis  of  data.  The  drawdown  curves  were  extrapolated  through  the  period  of  re¬ 
covery  and  the  recovery  measured  up  from  that  extrapolated  drawdown.  Values  of  these 
differences  in  h  ,  in  square  feet,  are  plotted  on  the  semi-log  graphs  in  Figures  D-39  through 
D-44.  These  graphs  are  os  follows: 

D-39  Semi-log  Drawdown  Graph  for  Aden  Well 


D-40 

Recovery 

Aden  Well 

D-41 

Drawdown 

PI,  P2and  P3 

D-42 

Recovery 

PI,  P2and  P3 

D-43 

Drawdown 

P4,  P5  and  P6 

D-44 

Recovery 

P4,  P5  and  P6 

Straight  lines  are  drown  through  the  points  plotted  on  the  semi-log  drawdown  and 
recovery  graphs.  From  the  slopes  of  these  straight  lines  the  hydraulic  conductivity  is  cal¬ 
culated,  and  from  their  intercepts  on  the  zero-drawdown  line,  the  storativity  is  estimated. 
Table  D-4  summarizes  the  data  from  this  interference  test  of  the  Aden  well.  The  first 

column  gives  the  name  of  the  observation  well  and  the  second  its  distance  from  the  pumping 

2  2 

well,  in  feet.  The  third  column  gives  the  change  over  one  log-cycle  of  h  ,  in  ft  . 

2  2 

The  average  change  of  h  over  one  log-cycle  is  11,1  ft  for  the  five  pairs  of  example 
of  drawdown  and  recovery.  Excluding  those  for  the  Aden  well  itself,  the  overage  is  9.9  ft'^. 
The  overage  intercept,  weighted  according  to  the  number  of  observation  wells  contributing 
to  each  number,  is  0,CX)26, 


TABLE  D-4 


DATA  FROM  INTERFERENCE  TEST  OF  ADEN  WELL 


[1] 

[2] 

[3] 

[4] 

Observat ion 
well 

Distance 

r/[ft] 

Change  over^^ 
one  log-cycle  of 
lA(h2)| 

1 ntercept 

[min/ft^] 

Aden 

_ 

15.0  d* 

_ * 

16.8  r* 

_ ♦ 

PI,  2,  3 

50.2,  100.2,  199.3 

13.0  d 

.0047 

8.4  r 

.0016 

P4 

51.3 

11.0  d 

.0035 

9.7  r 

.0024 

P5 

98.9 

8.5  d 

.0017 

10.0  r 

.0018 

P6 

200.6 

8.5  d 

.0017 

10.0  r 

.0009 

Average 

11.1 

- excluding  starred  numbers 

9.9 

.0026 

The  d  means  drawdown;  the  r  means  recovery. 

Average  initial  depth  of  flow,  11.2  ft 

Average  hydraulic  conductivity 

K  =  2.30  a/2irA(h^)  =  1,230  ft/day 

Average  initial  apparent  storativity 

S*  =  2.25  Kh*(t/r^)  =  2.25  x  1,230  ft/day  x  11.2  ft  x  .0026  min/ft^ 
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FIG.  D-42 
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3 

The  overage  discharge  was  173  gpm  or  33,300  ft  /day.  Thus  the  hydraulic  conductivity 
is  calculated  to  be  about  1,230  Ft/day. 

Interpretation^  Knowing  the  hydraulic  conductivity,  it  is  possible  to  substitute  back 
in  the  Formula  for  the  drawdown  and  calculate  the  effective  average  storativity  (S*),  The 
results  of  such  calculations  are  depicted  on  Figure  D-45,  which  is  a  graph  of  the  variation 
of  apparent  storativity  during  the  drawdown,  and  on  Figure  D-46,  which  is  a  similar  graph 
for  the  apparent  storativity  during  recovery.  Two  things  are  apparent:  first,  the  storativity 
increases  in  time  and  continues  to  increase  throughout  the  3-day  drawdown  period  and  5-day 
recovery  period;  second,  the  apparent  storativity  generally  is  greater  on  recovery  than  on 
drawdown.  Moreover,  the  apparent  storativity  is  greater  along  the  northeast  line  (Pi,  2  and 
3)  than  it  is  along  the  southeast  line  (P4,  5  and  6).  This  suggests  the  possibility  of  anisotropy 
as  regards  permeability,  for  we  arrived  at  these  storativities  by  using  a  uniform  hydraulic 
conductivity  calculated  from  the  average  change  over  one  log-cycle  of  drawdown  in  time 

(see  Column  [3],  Table  D-4).  The  slope  of  the  semi-log  time-drawdown  graph  for  PI,  2  and 

•  2  •  2 

3  is  13.0  ft  ,  while  that  on  the  semi-log  time  recovery  graph  is  only  8.4  ft  ,  the  average 

2 

being  10.7  ft  .  On  the  other  hand,  the  average  slope  of  the  semi-log  time-drawdown  graphs 

•  2  •  •  2 

for  P4,  5  and  6  is  9,3  ft  ,  and  of  the  semi-log  time-recovery  graphs  9,9  ft  .  As  the  difference 

in  apparent  storativity  on  the  two  lines  is  greater  than  the  difference  in  hydraulic  conductivity, 

it  is  evident  that  that  difference  in  apparent  storativity  is  not  due  to  our  assuming  a  uniform 

hydraulic  conductivity  for  both  lines.  Rather  it  would  appear  that  the  degree  of  confinement 

and  freedom  of  drainage  is  different  in  the  general  neighborhood  of  PI,  2  and  3  from  what  it 

is  in  the  vicinity  of  P4,  5  and  6 , 

Table  D-5  summarizes  the  variation  of  apparent  storativity  shown  graphically  on  Figures 
D-45  and  D-46.  After  4,000  minutes,  or  less  than  three  days,  the  apparent  storativity  on 
the  northeast  line  had  reached  about  17  percent  and  was  still  rising,  whereas  that  on  the 
southeast  line  had  reached  about  10,7  percent.  At  the  end  of  about  7,000  minutes,  or  less 
than  three  days  of  recovery,  the  apparent  storativity  on  the  northeast  line  had  exceeded  26 
percent  and  was  still  rising  sharply,  whereas  that  on  the  southeast  line  had  reached  something 
more  than  8  percent  and  was  rising  less  rapidly. 

It  is  now  evident  that  in  order  to  be  able  to  estimate  ultimate  storativity  or  "specific 
yield"  it  would  be  necessary  to  measure  interference  over  these  distances  for  a  few  weeks. 

Such  tests  are  not  now  feasible  but  may  become  so  during  the  decontamination  operation  in 
the  Future.  Pending  that  we  may  hazard  a  guess  that  on  the  northeast  line  the  ultimate 
apparent  storativity  would  reach  35  percent  and  perhaps  half  that  on  the  southeast  line. 


TABLE  D-5 


VARIATION  OF  APPARENT  STORATIVITY 
IN  TIME,  NEAR  ADEN  WELL 


[1] 

[2] 

[3] 

[4] 

[5] 

Storativity 

dur ina 

T  ime 

t 

Drawdown 

NLJ.me 

SE  1  ine 

Recove 

NE  line 

ry 

SE  line 

[min] 

(PI,  2,  3) 

(P4,  5,  6) 

(PI,  2,  3) 

(P4,  5,  6) 

100 

0.053 

— 

0.038 

0.052 

1,000 

.059 

.053 

.082 

.049 

2,000 

.092 

.067 

.121 

.058 

4,000 

.170 

.107 

.206 

.076 

7,000 

— 

— 

.266 

.083 

Est imated 

ultimate 

.18 

.35 

FIG.  D-46 
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FIG.  D-47a 


MAP  OF  PUMPING  TEST  LAYOUT  FOR  ADEN  WELL 
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Interference  Test  of  Marty  vVell 


Theory  and  procedure .  Same  as  before 

Conduct  of  test.  The  Marty  Well  was  pumped  continuously  for  about  4/4CX) 
minutes.  The  recovery  was  measured  at  frequent  intervals  for  another  3^300  minutes 
and  daily  thereafter  for  several  days. 

Three  piezometers  were  drilled  on  each  of  two  lines^  one  bearing  approximately 

southeast  and  the  other  one  approximately  southwest.  (See  map  on  Figure  D-48a.)  The 

distances  from  the  Marty  Well  to  the  piezometers  are  given  in  column  2  of  Table  D-6, 

Water-level  readings  made  in  the  Marty  Well  and  six  piezometers,  converted 

2  2 

to  read  in  terms  of  the  square  of  the  depth  of  flow  (h  ),  in  ft  ,  are  plotted  on  Figures 
D-49a  through  D-55b. 

Analysis  of  data.  The  drawdown  curves  were  extrapolated  through  the  period  of 

recovery  and  the  recovery  measured  up  from  that  extrapolated  drawdown.  Values  of 

2 

these  differences  in  h  ,  in  square  feet,  are  plotted  on  the  semi-log  graphs  in  Figures  D-56 


through  D-61 .  These  graphs  are 

as  follows: 

D-56  Semi 

-log  Drawdown  Graph  for 

Marty  Well 

D-57 

Recovery 

Marty  Well 

D-58 

Drawdown 

P3,  P2  and  PI 

D-59 

Recovery 

P3,  P2  and  PI 

D-60 

Drawdown 

P4,  P5  and  P6 

D-61 

Recovery 

P4,  P5  and  P6 

Straight  lines  are  drawn  through  the  points  plotted  on 

the  semi-log  drawdown 

and  recovery  graphs.  From  the  slopes  of  these  straight  lines  the  hydraulic  conductivity  is 
calculated,  and  from  their  intercepts  on  the  zero-drawdown  line,  the  storativity  is  es¬ 
timated. 

Table  D-6  summarizes  the  data  from  this  interference  test  of  the  Marty  Well. 

The  first  column  gives  the  name  of  the  observation  well  and  the  second  its  distance 
from  the  pumping  well,  in  feet.  The  third  column  gives  the  slope  of  the  line  on  the 
semi-log  graph,  in  feet  squared. 

The  overage  discharge  during  the  test  was  about  638  gpm  or  123,000  ft^/day. 
From  the  known  discharge  and  from  the  slopes  given  in  column  3,  values  of  hydraulic 
conductivity  in  column  4  were  obtained.  These  range  from  1,050  to  1,800  ft/doy  and 
average  about  1,330  ft/doy. 

The  fifth  column  of  Table  D-6  gives  the  intercepts  on  the  zero-drawdown  line  in 
min/ft  .  From  these  intercepts  values  of  the  ratio  h*/S*  in  feet,  given  in  column  6, 
were  obtained.  The  symbol  h*  stands  for  the  effective  average  depth  of  flow.  At  any 
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TABLE  D-6 

DATA  FROM  INTERFERENCE  TEST  OF  MARTY  WELL 


[1] 

[2] 

[3] 

[4] 

[5] 

[6] 

[7] 

[8] 

Observat ion 

Distance 

§lop®o  . 

Hydraul ic 
conductivity 
K/[ft/day] 

1 ntercept 
^(t/r^lo 

h*/S* 

h* 

s* 

wel  1 

r/[ft] 

A,hV[ft"] 

[min/ft^] 

FtT 

IfH 

Marty 

— 

30  r 

1,500 

— 

— 

P3 

64.5 

25  d 

1,800 

0.0016 

220 

13.8 

.063 

P2 

115.5 

34  d 

1,330 

.0018 

260 

19.5 

.075 

P1 

213.0 

29  d 

1,560 

.0017 

240 

22.5 

.094 

P4 

60.8 

31  d 

1,450 

.00081 

540 

17.5 

.033 

P5 

112.6 

35  d 

1 ,290 

.00100 

500 

20.7 

.041 

P6 

212.5 

31  d 

1,450 

.00125 

350 

21 .1 

.060 

P3 

64.5 

29  r 

1 ,560 

.0043 

95 

13.8 

.145 

P2 

115.5 

42  r 

1,070 

.0039 

150 

19.5 

.130 

PI 

213.0 

42  r 

1 ,070 

.0031 

190 

22.5 

.118 

P4 

60.8 

35  r 

1,290 

.0021 

240 

17.5 

.073 

P5 

112.6 

43  r 

1,050 

.0027 

225 

20.7 

.092 

P6 

212.5 

43  r 

1,050 

.0027 

225 

21.1 

.094 

SE  Line  d 

1,560 

240 

.077 

SW  Line  d 

1,400 

460 

.045 

SE  Line  r 

1,230 

145 

.131 

SW  Line  r 

r,130 

230 

.086 

Average 

1,330 

270 

.085 
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one  point  of  observation  this  cannot  exceed  the  initial  depth  of  flow.  To  determine  the 
limits  on  the  apparent  storativity,  the  initial  depths  of  flow  are  used  In  place  of  h*  in 
column  7, 

In  column  8  values  of  apparent  overage  storativity  are  calculated  from  the  figures 
given  in  column  6  and  7.  These  range  from  3,3  to  14,5  percent  and  average  about  8,5 
percent. 

At  the  bottom  of  Table  D-6  averages  are  given  for  each  of  the  two  lines  of  pie¬ 
zometers,  both  on  drawdown  and  recovery. 

It  should  be  noted  that  the  apparent  storativity  is  smaller  on  both  lines  for  draw¬ 
down  than  It  is  for  recovery. 

Interpretation.  The  storativity  remained  essentially  constant  throughout  the  period 
of  drawdown,  and  again  remained  essentially  constant  though  differing  in  value  throughout 
the  subsequent  period  of  recovery.  By  its  constancy  and  low  value  the  storativity  indicates 
that  there  may  be  fine-grain  sediments  in  the  upper  part  of  the  aquifer  overlying  the  more 
permeable  gravel.  In  which  there  is  capillary  retention  of  water.  An  examination  of  the 
logs  of  the  wells  verifies  that  this  condition  exists. 
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Pumping  Test  of  Nesom  Well  In  2S67W22 


Well  22bcc1  In  Section  22  of  T2S^  R67W  was  tested  by  the  U.  S.  Geological 

Survey  In  1955  or  1956.  This  well  Is  47.6  feet  deep  and  Is  cased  with  18-Inch  casing. 

It  was  drilled  In  1950  and  Is  owned  by  Harry  Nesom.  The  measuring  point  for  water 

levels  Is  0,8  feet  above  land  surface.  When  measured  13  September  1955, the  water 

level  stood  33,67  feet  below  the  measuring  point.  The  measured  yield  was  305  gpm 

with  a  drawdown  of  3.8  feet  after  10  hours.  (See  Petri  [1956],  Table  13.) 

By  the  Thiem  Method  the  transmissivity  was  found  to  be  about  128,000  gpd/ft, 

2 

or  about  17,000  ft  /day.  They  used  a  saturated  thickness  of  15  feet  and  got  a  hydrau¬ 
lic  conductivity  ("coefficient  of  permeability")  of  8,500  gpd/ft^.  We  calculate  14 
feet  for  the  saturated  thickness  and  get  about  1 ,300  ft/day  for  the  hydraulic  conductivity. 


Test  of  RMA  Well  In  3S67W4 


Well  4bca1  In  Section  4  of  T3S,  R67W  was  tested  by  the  Corps  of  Engineers  In 
January  1952  [U.  S,  Army,  1953,  p  3].  This  well  was  constructed  In  September  1938 
ard  was  cased  with  60-Inch  galvanized  casing  to  a  depth  of  51  feet  and  with  48-  Inch 
galvanized  casing  from  51  to  96  feet.  The  48-Inch  casing  was  gravel  packed  and  slotted. 

(Petri  [1956],  Table  13,  shows  24  Inches  for  a  diameter  of  this  well,  but  this  Is  evidently 
an  error,) 

The  well  was  tested  for  44  hours  and  produced  1 ,000  gpm  with  a  drawdown  of  22.3 
feet.  During  the  last  six  hours  the  pumping  rate  was  Increased  to  1 ,200  gpm  and  the  draw¬ 
down  was  32.  7  feet.  Observations  of  drawdown  were  mode  In  two  observation  wells^50 
feet  away.  In  one  200  feet  away,  and  In  one  300  feet  away,  but  the  data  were  Incomplete. 

The  well  was  tested  again  19  July  1953  for  64  hours  and  20  minutes.  The  discharge 
rate  was  500  gpm.  Observations  were  mode  of  the  drawdown  In  several  nearby  observation 

wells.  (See  USA  [1953]  Plate  1, )  The  data  were  analyzed  by  the  Thiem  Method,  The 

2 

average  transmissivity  was  found  to  be  136,000  gpd/ft  or  about  18,000  ft  /day.  The 
storativlty  was  found  to  be  .052.  The  averages  of  the  results  by  the  Thels  graphical  method 
were  130,000  gpd/ft,  or  17,000  ft^/day,  for  the -transmissivity,  and  .06  for  the  storativlty. 

As  the  Initial  saturated  thickness  was  45  feet,  the  hydraulic  conductivity  was  about  380  ft/day. 

Characteristics  of  Wells 


RMA  Well  In  3S  67W4.  Using  the  values  of  transmissivity  and  storativlty  determined 
from  the  Corps  of  Engineers  pumping  test  of  this  well,  and  taking  the  effective  radius  of  the 
well  to  be  24  Inches  or  2.0  feet,  the  drawdown  at  the  face  of  the  well  after  one  day's  con¬ 
tinuous  pumping  can  be  calculated  as  follows: 
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s  =  (2.30  Q/47TT)log(2.25Tt/Sr^^) 

=  0.93  ft  X  5.2  =  4.8  ft. 

That  is,  the  drawdown  at  the  face  of  the  well  is  about  4.8  feet.  This  may  be  compared 
with  the  drawdown  on  the  graph  on  Plate  8  of  the  report  by  the  U.  S,  Army  [1953].  The 
top  inclined  line  marked  "overage  observed  gradient"  intersects  the  vertical  line  for  a 
distance  of  2.0  feet  (representing  the  radius  of  the  well)  at  5,5  feet.  These  two  figures 
are  in  apparent  good  agreement,  when  allowance  is  made  for  the  difference  in  elapsed 
time.  However,  as  indicated  in  Appendix  B  under  the  subsection  entitled  "Tests  In  thin 
confined  aquifers,"  the  Theis  theory  must  be  modified  to  apply  to  h  rather  than  to  s. 
Moreover,  allowance  much  be  mode  for  well  losses  that  increase  with  drawdown. 

Discharge-drawdown  relation.  Figure  D-62  is  a  discharge-drawdown  graph  for 
Well  4bca1  at  Rocky  Mountain  Arsenal.  Recorded  thereon  are  reported  measurements 
of  discharge  and  drawdown.  In  January  1952  the  drawdown  was  measured  at  two  different 
rates  of  discharge.  During  the  pumping  test  of  July  1953  the  discharge  was  measured. 

The  formation  loss  at  that  rate  of  discharge  can  be  inferred  from  Plate  8  of  the  Corps  of 
Engineer's  report  [U.  S.  Army,  1953].  The  U,  S.  Geological  Survey  reported  a  yield 
of  650  gpm  at  a  (calculated)drawdown  of  about  21  feet  on  September  26,  1955  (see  Petri 
[1956,  Table  13]). 

The  bottom  of  the  formation  is  indicated  at  a  drawdown  of  45  feet.  The  scale  on 
the  right  hand  side  of  the  graph  gives  the  head,  or  the  depth  of  flow  in  feet,  measured 
above  the  bottom  of  the  formation. 

The  upper  curved  line  in  Figure  D-62  represents  the  calculated  formation  loss  at 
different  values  of  discharge  after  one  day  of  continuous  pumping.  The  lower  curved 
line  drawn  through  the  two  plotted  points  represents  an  estimate  of  the  discharge-draw- 
down  relationship.  It  is  seen  that  the  maximum  capacity  of  the  well  is  apparently  about 
1,350  gpm. 

The  curve  giving  the  combined  drawdown,  or  sum  of  formation  loss  and  well  loss, 
was  determined  by  assuming  that  the  well-loss  coefficient  is  inversely  proportional  to 
the  square  of  the  depth  of  flow.  The  point  obtained  on  Sept.  26,  1955  falls  some  dis¬ 
tance  off  the  curve.  The  discharge  was  only  reported,  not  measured.  This  may  account 
for  part  of  the  discrepancy.  The  rest  may  be  due  to  regional  drawdown  between  1952  and 
1955. 
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Symbols 


r  =  distance  from  pumping  well 
=  effective  radius  of  pumping  well 
t  =  time  elapsed  since  start  of  pumping 
h  =  depth  of  unconfined  flow  at  r  and  t 
h*  =  effective  overage  depth  of  flow 
Q  =  discharge  of  well 
K  =  hydraulic  conductivity  of  formation 
S*  =  effective  average  storativity 
C  =  well -I OSS  coefficient 


FIG.  D-62 
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DIFFERENTIAL  EQUATIONS  FOR  NONSTEADY 
UNCONFINED  GROUNDWATER  FLOWS  AND 
THEIR  DIFFERENCE  APPROXIMATIONS 

Homogeneous  Isotropic  Aquifer  in  Rectangular  Coordinates 

In  rectangular  coordinates  the  second-order  partial  differential  equation 
for  nonsteady  unconfined  flow  of  groundwater  in  a  homogeneous  isotropic  aquifer 
with  a  sloping  or  curved  bottom  is 

^  lx]  1^]  =  -K  ft 

(For  the  meaning  of  these  and  other  symbols  see  the  Glossary  at  the  end 
of  the  last  Appendix.) 

In  the  general  case  z  =  z  (  x,  y) .  That  is,  x  is  a  function  of  both  x  and  y  . 
Equation  [l]  may  be  reduced  to  universal  or  nondimensional  form  by  putting 
x-*4h*,  y -*  y  h^,  z ^  ho  ond  ht,  ,  also  i Sh^/K ,  Thus, 

Homogeneous  Isotropic  Aquifer  in  Cylindrical  Coordinates 
In  cylindrical  coordinates.  Equation  [l]  becomes 

Putting  t -*  and  transforming  the  other  symbols  as  before. 

If  the  bottom  of  the  aquifer  is  a  sloping  ptane,  we  may  write  z  or 

^  =  oT/O  c/Xf  6  ^  where  ct  is  the  slope  and  A  is  the  azimuth. 

If  the  bottom  of  the  aquifer  is  a  horizontal  plane  ,  Equation  [4j  may  be  sim¬ 
plified  to 


E.2 


PutHng  AiyO  =  X  this  becomes 


which  is  more  convenient  for  numerical  integration. 
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Inhomogeneous  Anisotropic  Aquifers 

If  the  aquifer  is  "orthotropic”  as  regards  hydraulic  conductivity,  then 
Kx  *  Ky.  If  it  is  also  inhomogeneous,  SKy,  /9x  ^  o  and  PKy /  o  , 

In  this  case,  therefore, 

■'“ly  ]  *  5fr 

By  putting  and  Ky  “♦  *^7  VK^  Ky  and  transforming  the 

other  symbols  as  before. 


Homogeneous  Anisotropic  Aquifers 

If  the  aquifer  is  orthotropic  but  homogeneous  as  regards  hydraulic  conductivity, 
then  =o  and  =  o  .  Putting  I and  ”7 , 

Equation  [8]  becomes 

[(!(- =  ;^  [9j 

This  transformation  preserves  the  area:  L7)  -  A|  A^,, 


Difference  -  Differential  Equations 

The  derivation  of  the  difference-differential  equations  for  use  in  "passive  " 
RC  networks  (direct  analogs)  or  in  DC  analog  computers  (indirect  analogs)  is  best 
seen  by  the  special  case  of  Equation  [2],  restricted  to  one  space-coordinate,  with 


^  i'i^)  ^  -L 

?|  1  dr 


[loj 


We  use  i  as  the  index  of  nodes  in  the  mesh,  with  spacing  A  f  .  Thus 
Equation  [10]  becomes 

If  there  are  *n  nodes,  there  are  m  simultaneous  non-linear  difference- 
differential  equations  like  [11]  to  be  solved.  This  the  analog  computer  or  the 
passive  network  does  automatically  if  the  modeling  is  appropriately  carried  out. 


Difference  Equations 


For  use  with  a  desk  calculator  or  electronic  digital  computer  the  equations 
are  differenced  both  in  time  and  in  space.  We  may  use  j  as  the  index  of  points 
in  time  and  employ  the  notation  ~  ~ 

using  subscripts  for  space  and  superscripts  for  time,  understanding  that  for 
simplicity  either  index  may  be  omitted.  Thus, 

iii-i  -  %( j)- 

V*'  = 


With  this  notation  Equation  [10]  becomes 

AV'' =  ^<*1  ~  ^  ^  <-/  _ 

by  which  the  m  values  of  can  be  calculated  from  the  rn  known  values  of 

Similar  difference  approximations  may  be  set  up  for  Equations  [2],  [4],  [6], 
[8]  and  [9]. 
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UNCONFINED  RADIAL  FLOW  TO  WELLS 
Single  Well  in  Effectively  Infinite  Aquifer 

The  non-linear  partial  differential  equation  for  nonsteady  unconfined 
flow  to  a  single  isolated  well  has  already  been  given  (Equation  E-[4j).  If  the 
aquifer  is  effectively  infinite,  that  "boundary  condition"  may  be  expressed  as 
follows:  {9%  ^ ^ 

Boundary  Conditions 

As  the  differential  equation  is  of  second  order  \n  ^ ,  a  second  (inner) 
boundary  condition  in  f  must  be  specified  to  completely  define  the  problem. 

For  example,  a  "step"  well  with  full  instantaneous  drawdown  would  have 
=  o  for  for  all  positive  ^  .  A  "ramfi' well  might  attain 

full  drawdown  at  a  uniform  rate  during  a  short  initial  period  of  time.  A  well 
of  constant  discharge  would  have  O  /P\)^  »  const,  for  all  positive  T  . 

Initial  Condition 

As  the  partial  differential  equation  is  of  first  order  in  T,  one  condition 
on  ^  in  T  (usually  an  initial  condition)  must  be  specified.  In  the  foregoing 
examples  the  initial  condition  is  ^  (X)  =  o  atT  ^  O  for  all  X  . 

Bounded  Aquifer 

If  there  is  a  single  well  in  a  bounded  horizontal  aquifer  with  no  inflow 
from  outside  or  a  large  number  of  wells  uniformly  spaced  on  a  regular  pattern 
In  such  an  aquifer  the  outer  boundary  conditions  may  be  expressed  as  ^  / Pv  - 

on  the  outer  boundary,  where  V  is  the  outward  directed  normal. 

Wells  in  Square  Array 

An  approximate  solution  for  the  nonsteady  flow  to  wells  in  uniform  square 
array  may  be  obtained  by  replacing  each  square  tributary  to  a  single  well  by  a 


circle  having  the  same  area,  whose  radius  therefore  is  l^e  =  2a/-\/F  ,  2a  being 
the  spacing  of  the  wells.  An  exact  solution  for  steady  flow  may  be  obtained  in 
terms  of  the  Green's  function  for  a  square.  An  "exact"  solution  for  nonsteady 
flow  may  be  obtained  only  by  analog  or  fine-mesh  digital  computation. 


Wells  in  Hexagonal  Array 

Similarly,  approximate  solutions  for  the  nonsteady  flow  to  wells  in  uniform 
hexagonal  array  may  be  obtained  by  replacing  each  hexagon  tributary  to  a  single 
well  by  a  circle  of  the  same  area. 


Wells  Nonuniformly  Spaced 

The  nonsteady  flow  to  wells  unevenly  spaced  in  an  unconfined  aquifer  with 
a  sloping  or  curved  bottom  is  best  handled  either  by  a  passive  network  or  by  an 
electronic  analog  computer.  Equation  E-[2]  may  be  used,  with  appropriate  modi¬ 
fication  for  locally  radial  flows  and  drawdowns  at  the  wells. 


Line  of  Wells 


The  steady  flow  of  water  to  a  typical  well  of  discharge  Q  and  radius  in  an 
infinite  line  of  uniformly  spaced  wells  is  given  by  the  equation 

-fr^n/4a.^ 


^  ■  irrK  ^ 


where  h^  is  the  head  at  the  well-face  and  2  a.  is  the  well  spacing. 

An  approximate  solution  for  the  nonsteady  flow  to  a  line  of  wells  may  be 
obtained  in  terms  of  the  Green's  function  for  an  infinite  strip  with  impermeable 
boundaries.  The  "exact"  solution  may  be  obtained  only  by  fine-mesh  analog  or 
digital  computation. 

Line  of  wells  rrf>«- grad  lent.  If  the  line  of  wells  runs  cross-gradient  and 
if  the  initial  depth  of  (inclined)  flow  is  uniform,  the  flow  may  be  described  by 
solving  Equation  E-[2]  with  H,  -  cC  and  with  the  conditions  on  the  "impermeable" 
boundaries  or  dividing  stream-lines  x  =  ±  Q.  appropriately  specified:  Qh/3x.  =  0  . 
In  this  case  Equation  E-[2]  is  solved  by  fine-mesh  analog  or  digital  computation. 

Line  of  wells  down-gradient.  If  the  line  of  wells  runs  down-gradient,  even 
if  the  initial  depth  of  (inclined)  flow  is  uniform,  still  the  flow  cannot  be  divided 


by  parallel  lines  as  before.  In  this  case  the  solution  is  found  for  the  infinite 
plane  by  fine-mesh  analog  or  digital  computation. 

Wells  In  Square  Array  in  inclined  Uniform  Aquifer 

With  the  wells  uniformly  spaced  in  square  array  in  an  aquifer  of  initially 
uniform  thickness  or  depth  of  flow  the  flow  can  be  divided  by  parallel  streamlines 
running  down-gradient.  Again  Equation  E-[2]  is  solved  by  fine-mesh  analog  or 
digital  computation  with  ^^  =  of>7  and  with  (9h/Px)  =  0  af  y  =  ±  a  . 


UNCONFINED  UNIDIRECTIONAL  FLOW  TO  TRENCHES 


Horizontal  Aquifer 


The  differential  equation  for  the  nonsteady  flow  of  groundwater  to  a 
trench  cutting  through  a  thin  horizontal  unconfined  aquifer  has  already  been 
given  (Equation  E-[2] ).  With  the  trench  running  in  the  ">7 -direction  the 
equation  becomes 


£11!  -  £l 

9'r 


The  corresponding  difference  equation  is 


z 


2  A|*/AT 


2 


For  stability  we  choose  AT  /  2  . 

Figure  G-1  is  a  logarithmic  graph  of  the  water-table  decline  near  a 
trench  with  full  drawdown  in  an  aquifer  with  horizontal  bottom  and  initially 
horizontal  water-table.  The  plotted  points  were  obtained  by  numerically 
integrating  Equation  [2]  using  a  desk  calculator.  Only  a  limited  range  of  the 
composite  variable  |/2  Vt  is  covered.  In  order  to  calculate  the  variation  of 
flow  into  the  trench,  a  much  longer  run  on  a  digital  computer  with  finer  mesh 


is  needed. 


Inclined  Aquifer 

The  differential  equation  for  the  nonsteady  flow  of  groundwater  to  a  trench 
cutting  through  a  thin  inclined  unconfined  aquifer  has  also  been  given  (Equation 
E-[2j,  with  -  4  ).  Putting  of* T  -;»T'weget 

-||7  [( "^  -  I 

Water-table  decline  upstream.  Figure  G-2  is  a  logarithmic  graph  of  the 
water-table  decline  upstream  from  a  trench  with  full  drawdown  in  an  aquifer 
with  inclined  plane  bottom  and  initially  parallel  inclined  water  table.  The 


plotted  points  for  different  values  of  4  were  obtained  by  numerically  integrating 
Equation  [3]  using  a  desk  calculator.  Again,  only  a  limited  range  of  the  composite 
variable  4  /2*^'  is  covered.  Digital  calculations  with  electronic  computer  are 
needed  to  cover  the  full  range  and  accurately  calculate  the  inflow  into  the  trench 
from  upstream. 

Water-table  decline  downstream.  Figure  G-3  is  a  similar  logarithmic  graph 
of  the  water-table  decline  downstream  from  such  a  -trench.  The  decline  is  much 
more  pronounced  downstream  than  upstream. 

Figures  G-1,  2  and  3  are  merely  Illustrative.  Nearly  instantaneous  draw¬ 
down  was  applied  at  the  face  of  the  trench.  Other,  more  realistic  examples  must 
be  run  to  elucidate  the  behavior  of  such  a  non-linear  system. 
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SEMI-CONFINED  FLOWS  AND  MIXED  FLOWS 


Transmissive  Aquifer  with  Storative  Semi-Confining  Layer 

Assume  a  relatively  transmissive  aquifer  of  thickness  Is  and  hydraulic  con¬ 
ductivity  K  is  overlain  by  a  much  less  transmissive  but  more  storative  confining 
bed  of  hydraulic  conductivity  K  and  storativity  S^.  The  flow  may  be  con¬ 
sidered  to  be  completely  vertical  in  the  confining  bed  and  completely  horizontal 
in  the  aquifer,  the  refraction  being  virtually  ?(?*at  their  interface. 

The  vertical  flow  in  the  confining  bed  is  governed  by  the  equation 


K'  ^ 


o'  9h 

^  0i 


where  h  is  the  head  in  the  aquifer  and  h^  is  the  elevation  of  the  water  table  in 
the  confining  bed. 

The  horizontal  flow  in  the  aquifer  is  governed  by 


JLk  ^  (  r  «  _  k'  H 

r  3r  ^  9rJ  ^  U 


9r  9ri  •'  y,' 

Equations  [l]  and  [2]  share  a  common  term,  which  relates  to  the  interflow 
between  the  confining  layer  and  the  aquifer. 

By  an  appropriate  reduction  to  nondimensional  form  Equations  [1]  and  [2j 
become 

I*' 

with  Xsj2n^. 

Equations  [3]  are  amenable  to  solution  by  digital  computation,  solving  the 
right-hand  equation  first  and  the  left-hand  equation  afterward  in  successive  alter¬ 
nations.  Their  equivalents  in  differenced  form  are: 

-  e-**-'  .7n+>).  ) 

For  a  steady  well  in  an  effectively  infinite  aquifer  the  inner  boundary  con¬ 
ditions  are  ^)o  *  const. ,  0  and  the  outer  boundary 

cond i tons  are  *  0=  .  The  initial  condition  might  be 

*  const. 


H-2 


Mixed  Flows 

Of  great  practical  interest  in  the  present  investigation  is  the  case 
where  a  radial  flow  may  initially  be  confined  or  semi-confined  uniformly 
over  all  p  and  where  unconfinement  may  develop  at  the  well  and  spread 
outward  from  the  well .  From  the  we  1 1 -face  (oJt  r  -  r^)  out  to  some 
moving  cylindrical  boundary  (at  r  «  r^)  the  flow  may  be  governed  by 
Equation  E-[6]  while  beyond  that  moving  boundary  it  may  be  governed  by 
Equation  [3], 

Crude  attempts  hove  been  mode  to  solve  this  system  of  equations  with 
little  success.  Computation  by  analog  or  digital  methods  will  give  satisfactory 
results. 


DISPERSION  OF  CONTAMINANTS  IN  GROUNDWATER 


Unidirectional  Dispersion 


Dispersion  in  a  uniform  unidirectional  flow  of  groundwater  is  governed 
by  the  equation 


_L 

V  9i 


where  0  is  the  dispersion  constant,  c  is  the  concentration,  and  V  is  the  overage 
velocity  (=  Q  /  Af,  with  -f  =  porosity) . 

The  difference  equation  corresponding  to  Equation  [Ij  is 


C  i-f  I  ~  2Ci  +  Ci. 


_  C  If  I  -  C  i-i 

2  Ax 


I  ACi 


This  is  amenable  to  analysis  and  has  a  solution  in  terms  of  the  error  function. 


Radial  Dispersion 

Dispersion  in  a  steady  radial  flow  of  strength  Q  in  an  aquifer  of  thickness  b 
is  governed  by 

n  3*c  _  9c  _  2-rrryf  9c 

dr  ~  Q 

This  equation,  or  rather  its  differenced  equivalent,  is  not  amenable  to 
exact  mathematical  analysis  but  should  be  solved  by  analog  or  digital  compu¬ 
tation. 


Stratified  Groundwater  Flow 

It  is  important  to  understand  the  theory  of  stratified  flow  in  unconfined 
groundwater  systems.  Consider  the  flow  of  fresh  water  of  viscosity /u  and  specific 
weight 'v  beneath  a  water  table  whose  elevation  is  h  in  an  aquifer  of  porosity  T 
and  permeability  k  .  Let  this  flow  be  bounded  beneath  by  its  interface  with  a 
parallel  or  opposing  flow  of  salt  water  of  viscosity//'  and  specific  weight  y' , 

The  elevation  of  the  Interface  is  z,  and  the  elevation  of  the  bottom  of  the  aquifer 
is  2'. 
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By  equating  pressures  at  the  interface  at  elevation  z  it  is  found  that  for 
"gradually  varied"  flow  y'h'  *  yh  +6^'-  y)^^  where  is  the  head  (assumed 
uniform  over  the  vertical)  in  the  salt  water. 

The  fresh-water  flow  is  governed  by  the  equation 


Compare  this  equation  with  Equation  E-[lJ. 

The  salt-water  flow  is  governed  by 

Using  the  relation  between  h  and  b  and  z  and  simplifying/  Equations 
[4]  and  [5]  may  be  written 


jS.  Sb.  —  ^  U  z )  .^ill  =  J-  Sjl  =  jfL.  ^  \  (z  -z'i  -S.  / /)  -  y  ~  z  )  1 
By  appropriate  reduction  to  nondimensional  form  Equations  [6]  become 


[6J 


3 


=  *  /I 


PT 


^■h  [7] 


2I 

9  T 

Equations  [7]  are  amenable  to  digital  computation  by  solving  the  left- 
harKl  equation  first  and  the  right-hand  equation  afterward/  in  successive  alter¬ 
nations.  Or/  they  may  be  solved  by  electronic  analog. 


Dispersion  in  Stratified  Flows 

The  dispersion  of  ions  in  stratified-flow  systems  may  be  described  by 
appropriate eombinations  of  the  concepts  formalized  under  the  three  preceding 
topics.  Also  allowance  may  be  mode  for  the  Influence  of  inhomogeneity  and 
anisotropy  by  incorporating  the  principles  of  Equations  E-[8]  and  E-[9]. 


J 


GLOSSARY 


Listed  and  explained  below  are  the  more  important  mathematical  symbols 
employed  in  the  appendices,  particularly  in  Appendices  E  through  L  For  other 
symbols  see  page  D-21 . 

O'  =  half-width  of  flow-strip  in  horizontal  plane, 
or  other  characteristic  horizontal  length. 

b  =  thickness  of  an  aquifer. 

C  =  concentration  of  an  ion  subject  to  dispersion  . 

D  =  dispersion  constant  of  a  porous  medium  such 
as  an  aquifer. 

f  =  porosity  of  aquifer . 

h  =  head,  or  elevation  of  water  table  above  arbitrary 
datum  (datum  is  usually  bottom  of  aquifer  in  hori¬ 
zontal  case).  Also,  depth  of  unconfined  flow. 

ho  =  head  at  x  =  oor  at  r  =  o  (origin  of  space  coordinates)  . 

h*  =  initial  head,  i.e.  at  t  =  o  (origin  of  time  coordinates). 

I  =  index  for  nodes  in  space-mesh  or  network. 

j  =  index  for  point  in  time  in  difference  equation. 

=  hydraulic  conductivity,  the  product  of  permeability 
by  specific  weight  of  fluid,  divided  by  viscosity  of 
fluid.  Commonly  miscalled  "permeability"  and 
"coefficient  of  permeability.' 

Kx  =  X- component  of  hydraulic-conductivity  tensor. 

=  y- component  of  hydraulic-conductivity  tensor. 

=  permeability 

Q  =  discharge  of  well 

S  =  storativity  or  "coefficient  of  storage"  (Theis) 

t  =  time  (usually  time  elapsed  from  start  of  pumping  or 
start  of  draining). 

V  =  average  velocity  of  groundwater  flow  (-  Q/^f) 

X  =  horizontal  coordinate,  normal  to  y 

y  =  horizontal  coordinate,  normal  to  X 

2  =  elevation  of  bottom  of  aquifer  above  arbitrary  datum,  z  =  ^  y ) 
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^  -  angle  or  slope  of  plane  bottom  of  aquifer. 

Tr  =  specific  weight  (weight  per  unit  volume)  of 
fluid  such  as  water  . 

=  nondimensional  vertical  space  coordinate,  or  head-ratio. 
'*}  =  nondimensional  horizontal  space  coordinate,  normal  to  4* 

ory/a. 

6  =  azimuthal  angle. 

^  -component  of  conductivity-ratio  tensor. 

Kyf  =  -component  of  conductivity-ratio  tensor. 

X  =  (natural  logarithm) 

=  viscosity 

P  =  nondimensional  coordinate,  normal  to  boundary 
^  =  nondimensional  horizontal  space  coordinate,  normal 

to  ■»7  .  I  =  >/ho  or  x/pL  . 

P  =  nondimensional  radial  coordinate  ^  . 

T  =  nondimensional  time  coordinate,  or  time-ratio. 

T  =  Kt/Shft  or  K  hot /Sa*^  depending  on  the  case. 
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NOTES: 

1.  FOR  RANGE  200-400  Mg/l  CHLORIDES  IN  CONTAMINATED  SOURCE  USE  CHART  1 

2.  FOR  RANGE  400-2200  Mg/L  CHLORIDES  iN  CONTAMINATED  SOURCE  USB  CHART  2 
a  FOR  RANGE  2200-5200  Mg/L  CHLORIDES  IN  CONTAMINATED  SOURCE  USB  CHART  5 

4.  THE  CHART  ISJASEO  ON  ZERO  CONCENTRATION  OF  CHLORIDE  IN  THE  RECEIVING  S  REAM 
^R.A<:rjAL  VALVE  OF  PERMISSIBLE  D.SCHARG5  FROM  CONTAMINATED  SOURCE  hTcf 'M  MULTIPLY  THE 
CHART  VALUE  &Y  THE  FACTOR  FOR  THE  MONTH  AND  SIREAM^OWN  IN  THE  TABLE  C'N_Tri IS  CHART. 

THE  FACTOR  IS  DETERMINED  FROM  THE  AVERAGE  CONCENTRATION  OF  CHLORIDE  CUrvINO  EACH  MONTH. 

5. FOR  EXAMPLE  Or  USE  OF  CHART,  SEE  CHART  1,  Fig.  lU 
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:ontaminated  source  use  chart  i 
contaminated  source  use  chart  2 

N  CONTAMINATED  SOURCE  USE  CHART  5 

ATION  OF  CHIORIDE  IN  THE  RECEIVING  STREAM 
iARGE  FROM  CONTAMINATED  SOURCE  O^M  V.ULTIP'Jf  THE 
DNTH  AND  STREAM  SHOWN  IN  THE  TABLE  ON  THIS  CHA"Rt. 
6RAGE  concentration  OF  CHLORIDE  CUHINO  EACH  MONTH, 
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NOTES:  ^ - 

1.  CLASSIFICATIONS  AND  DESCRIPTIONS  AT  THE  RIGHT  Cf  THE  BORING 
PROFILE  ARE  SUPPLEMENTAL  DATA  BASED  ON  FIELD  INSPECTION  OF  THE 
SAMPLES  UNLESS  OTHERWISE  NOTED. 

2.  ALL  elevations  REFER  TO  MEAN  SEA  LEVEL,  1929  GENERAl  ADJUSTMENT. 

3.  SOIL  CLASSIFIED  BY  THE  UNIFIED  SOU  CLASSIFICATION  SYSTEM. 


LEGEND: 

90  DRILL  HOLE  NUMBER 

ELEV.  5108.81  ELEVATION  TOP  OF  HOLE 

DATE:  MAY  60  DATE  ORUtlNG  COMPLETED 

•  INDICATES  LABORATORY  CLASSIFICATION 
CL  LEAN  OR  SANDY  CLAY 
CH  FAT  OR  SANDY  FAT  CLAY 
SP  SAND.  POORLY  GRADED 
SW  SAND.  WELL  GRADED 
CC  GRAVEL,  CLAYEY  this  0 

GM  GRAVEL.  SILTY 
SC  SAND.  CLAYEY 
SM  SAND.  SILTY 
GP  gravel,  POORLY  GRADED 
GW  GRAVEL,  WELL  GRADED 
Ml  SILT.  LOW  PLASTICITV 
MH  SILT,  HIGH  plasticity 
LI  LIQUID  LIMIT 
PI  PLASTICITY  INDEX 
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Piping  ond  well  ioyout  is  scHnmatic  ond  dons  nnt  nnenssority 
•xoci  position  of  '•'otIs  or  location  ct  pipot. 

Pipos  or*  sizod  to  occommedoto  Nioximwm  onticipetod  flo^*  os 


( 


WELL  flELD  A  2000  gpin 
WELL  FIELD  B  2000  gpm 
WELL  FIELD  C  1000  gpm 
SUMP  A  2000  gpm 
SUMP  »1  1000  gpm 

SUMP  B7  tOOO  gpm 
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Recovery  Tests 


MAP  OF  PART  OF  CONTAMINATED  / 
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D  AREA  NEAR  ROCKY  MOUNTAIN 
.RIZING  WELL  TESTS 
erence  Tests 


TABLE  D-1 


SUMMARY  OF  DRAWDOWN-RECOVERY  TESTS  OF  WELL! 


[1] 

[2] 

[3] 

1 - 1 

1 - 1 

[5] 

[6] 

[7] 

[8] 

[9] 

[10] 

Depth/^ 

feet 

Casing 

Pump 

Reported 

Measured 

Well 

Locat ion 

d  iam. 

D  iam. 

Horse 

discharge, 

gpm 

d i scharae . 

name 

2S  67W 

inches 

Kind 

i nches 

power 

gpm 

lO^ft/day 

Aden 

9  daa 

24.1 

48 

concr . 

- 

10 

400 

194 

37.4 

Dobbs 

9  ddb 

45 

48/24 

steel 

5 

— 

— 

42 

8.1 

Marty 

21  add 

52 

48 

concr. 

8 

20 

— 

765^ 

147. 0^ 

Masunaga 

10  dec 

39 

48- 

concr . 

6 

5 

300 

176 

33.9 

Mat sumo to 

16  ddc 

36.9 

48 

concr . 

6 

5 

400 

336 

64.7 

Mat sumo to 

16  ddd2 

36.94 

48 

concr . 

6 

5 

500 

384^“ 

74. 0^ 

Monson 

1 5  badd 

40 

18 

steel 

6 

10 

400 

386 

74.4 

Myers 

10  bded 

39.5 

48 

— 

3 

3 

150 

46 

8.9 

Powers 

1 5  ced 

40 

48. 

concr  • 

6 

10 

800 

264 

50.9 

Wolpert 

4  adb 

18.1 

36 

concr . 

— 

— 

250 

275 

53.0 

Average 

37.2 

287 

55.2 

- exc 1 ud i ng  starred 

items 

215 

41 .4 

a) 

Depths  shown  to  nearest  foot  are  reported  depths  of  drilling; 
depths  shown  to  tenths  or  hundredths  are  sounded  depths. 
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[10] 

red 

arae  _ 

^  f t^day 

[11] 

Drawdown  over  1 
log-cycle  in  time 
A,(h")/[n^] 

[12] 

Hydraul ic 
conduct i V i ty 
K/[ft/day] 

[13] 

Drawdown 
at  ICO  min 
A(h")/[ftn 

[14] 

Khll 

Afh^) 

[15] 

[16] 
Initial 
depth  of 
flow, 

hV[rt] 

37.4 

10.3 

1,330 

86.1  -  29.0 

5.55 

0.77  -  4 

9.3 

8.1 

2.5 

1 ,190 

131.8  -  113.5 

7.33 

.94  -  6 

11.5 

147.0^ 

18. 

3,000* 

377  -  76 

16.7 

.97  -  15* 

19.4* 

33.9 

12. 

1,040 

102.8  -  39.8 

5.25 

.96  -  4 

10.1 

64.7 

9.0 

2,640 

131 .8  -  44.0 

9.76 

.85  -  8 

11.5 

74.0* 

10. 

2,720* 

201 .9  -  54.0 

14.8 

.83  -  13* 

14.2* 

74.4 

12.5 

2,180 

148.9  -  76.0 

5.83 

.70  -  4 

12.2 

8.9 

2.5- 

1 ,3004 

12.25  -r  .05 

4.94 

.41  -  3 

3.5 

50.9 

18. 

1,040 

118 

6.55 

.66  -  5 

11 .8 

53.0 

17. 

1 ,140 

150  -  18 

7.75 

.50  -  6 

12.3 

55.2 

1,758 

.96  -  7 

11 ,6 

41 .4 

1,483 

.72  -  5 

10.3 

013-M2-601116  5/15  CEJ/epf 


013-M2-601116  6/15  CEJ/epf 


013-M2-601116  8/15  CEJ/epf 


013-W2-601116  9/15  CEJ/epf 


013-M2-601116  10/15  CEJ/epf 


>CE\JFFEl.  ^  ESScR 


F.NCRAVIMG  ,'5.54-3,  10  'i  10  TO  rMF  MALF  U 
,'iiEN  ORDr-iftifio  diA7s  c:.M  o« ,  on-vv/o.-c  fin  inAciN*;  v* 


j  ' ;  i 

1  :  ■ 

'  1  •  ;  ‘  j 

i  ^  i  1  1  ^  ^  ^  1  1  ■ 

1  :  i-  •  '  i 

'  i 

gak)  L 

SOL  i 

j  iitoQQ,.L.;,i.  i  itS^i  i  iSab-i  ■-'L.J.g6(i)o.J  i 

i_fiQQQj 

kEUFFF.L  &•  fISStR  CO. 
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.^34*3  io  'X  10  10  THE  H-Mf  !NOM 

s'Ait  ?ci.oft,  eT?Al^••(;C  L-n  TrAftiNa  f>*.r'e* 


v*HrM  nr\t>rRmc  stat^  ooloh,  ynw  •  c  traciwu  ;-^pgj 


^  i  I  1  uebo  ‘ 


y-4Yfc> 


ENGRAVING  354-3,  10  X  10  TO  THE  HALE  INCH. 

CHrN  5t<'JErMNC  STATE  COLOR ,  OnAWINO  Oft  T.ftACiNO  PApEft. 


F!G.  D-50a 


RNGPAVING  10  X  10  TO  THt  HAtF  INCH 


D-50b 


ENGRAVING  33A  !0  X  10  TO  THE  HALF  INCH, 

/HEN  OROtRI/IG  STATE  CO  LOR ,  PRAWiNG  OR  TPACtHO  PAPFR 


,=',G.  D-51b 


KcUrPfcL  &  ESSER  CC  N  y. 


=\G.*D-52a 


ti:  RNIMVWO  '  UOTOO  3i»\..S  N-JH/ 

HONl  IH.I  O.l  01  X  Cl  t  otjiv/yon^ 


FEftENCE  OF ‘MARtYiWtLL 


Co^td . 


NOMI  JIVH  HHi  Oi  01  X  0  ONlAVtIONJ 


ENGRAVING  3S4-3,  10  y  10  TO  THE  HALF  INCH 

WHEN  ORO£«IMO  STATE  UOLOR,  DRAWING  OR  TflAOIKQ  v>Af>E! 


ENGRAVING  33/1  10  X  10  10  TI-!F  HAl  F  INCH 


KEUFFEt  &  ESSER  CO.  Y. 


EHGRAViNG  334-3,  X  10  TO  THE  KUr  INCH. 

V/HCN  ORf>fH!NC  STATE  CO  LOP ,  OHAWtMG  Oft  TftAClMO  kAPER 


FIG.  D-54b 


KEUFPeL  &ESSER  CO.  N  Y. 


feNGRAVlNG  334-3,  .J  X  10  TO  THE  HALF  INCH. 

WHEN  OKDCniNrt  STATE  OOt.Oft ,  ORAWIMC  Oft  TftAOIHO  PAPER. 


FIG,  D-55a 


rMGRfVING  .^34-3  X  10  TO  THK  HALF  INCH. 
(jpar-.Rific  htme  color,  on  thaojho  PAf-ec 


KEUFFEL  &  ESSER  CO.  N  Y. 


GRAPH  OF  WATER-TABL 
FULL  DRAWDOWN 

=30  o 


\BIE  DECLINE  NEAR  TRENCH  WITH 


N  IN  HORIZONTAL  AQUIFER 
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GRAPH  OF  WATER  TABLE  DECL 
WITH  FULL  DRAWDOWNLI 


gSamiSggiiiiisyHai 


GRAPH  OF  WATER  TABLE  DECLI 
WITH  FULL  DRAWDOWN 


CLINE  DOWNSTREAM  FROM  TRENCH 


TYPICAL  OBSERVATION  WELL  INSTALLATION 
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ce  USE  CHAft7  1 
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>u9ce  USE  CH-^sr  3 


HE  deceiving  S7S5A.V.. 

JEO  SOURCE  IN  Gr^  .•ULTIPLt  Tr-E 
*4  IN  THE  TAiilE  QN  ,5  CHART. 

OF  CMlORiCc  DURING  EACH 


FACTOR 

MOMTH  PLATTE  j  IRRIGATION  I 
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